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Preface 


The purpose of this book is to provide an introduction to analytic number 
theory. The theme is the mean value of certain important elementary arithmetic 
functions, the Euler <j> function, the divisor function, the squarefree function, 
and the prime divisor function. 

Analytic number theory, as the name indicates, is a blend of analysis and num¬ 
ber theory. It soon transpires that one needs algebra and geometry. 

Let us sketch the contents of the various chapters. In Chapter 1, various esti¬ 
mates used throughout the volume are given. In no other field of mathematics is 
it as important to have precise estimates. In Chapter 2, various transform 
methods are given, Fourier series both infinite and finite, the Laplace transform 
and the Mellin transform. The Mellin transform plays a very important role in 
number theory, as is indicated. In Chapter 3, some properties of congruences 
are given. In Chapter 4, some properties of the gamma function are given. In 
Chapter 5, the Riemann zeta function is introduced. This function plays many 
roles in number theory. It serves as a generating function for many arithmetical 
functions of interest and plays an important role in estimating error terms. In 
Chapter 6, the Poisson summation formula is discussed. This is one of the most 
powerful methods of obtaining important identities. In Chapter 7, the topic of 
functional equations is introduced. The utility of this method is shown by appli¬ 
cation to the theta functions. 

In Chapter 8, we turn to various mean values associated with Euler 0 function. 
As we show, a representation for this function makes the mean values easy to 
obtain. In Chapter 9, we consider analogous questions for the divisor function. 
Here, the corresponding questions are much more difficult. 

In Chapter 10, we consider some problems associated with squarefree numbers. 
What is surprising is how difficult these questions appear to be. 

Chapter 11 consists of three unrelated topics. In the first part of this chapter, 
we consider the prime divisor function, assuming the prime number theorem. In 
the second part of the chapter, we present the sieve method of Selberg. Finally, 
in the third part of the chapter, we present some results on the algebraic inde¬ 
pendence of the elementary arithmetic functions. 


xv 



XVI 


Preface 


Finally, in Chapter 12, we show the utility of Tauberian theory in analytic 
number theory. 

Much has been omitted. We have done this for several reasons. In the first 
place, there is no point to repeating material which is readily available. In the 
second place, it would not be in anybody’s interest to make the volume too large. 
In the third place, we have concentrated on the areas in which we have had 
greatest personal interest. The connections with algebra, particularly, algebraic 
number theory and group theory, have been hinted at. Many other connections 
are referred to—either in the text or bibliography and comments. 

Any serious student of number theory will want to refer to the following 
books: 

Dickson, L. E., History of the Theory of Numbers, Chelsea, New York, 1952. 
Landau, E., Vorlesungen iiber Zahlentheorie, Teubner, 1921. 

Hardy, G. H., and E. M. Wright, An Introduction to the Theory of Numbers, 

Clarendon Press, Oxford, 5th ed., 1979. 

It will be clear how much has been omitted here. In particular, we do not dis¬ 
cuss the important subject of trigonometric sums. This topic has been omitted 
for two reasons. In the first place, the discussion is very technical which makes it 
unsuitable for an introductory text. In the second place, no brief account of this 
important area can be given. We have noted a reference to an excellent exposi¬ 
tory paper by Vinogradov in which many other references will be found. 

It will be evident from the text how many open questions there are in analytic 
number theory. It is thus an ideal field for research. 

We presuppose a knowledge of elementary number theory and the rudiments 
of the theory of functions of a complex variable. Some of the miscellaneous 
exercises require a bit more. With this background, the book can be used for a 
one-semester course, or self-study. 
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Analytic Number Theory 

An Introduction 


1. Estimates 


1. Introduction 

In this chapter, we will obtain some results which we shall use in the rest of 
the book. 

First we shall compare the sum in Xff =0 f(k) with the integral f%f(x)dx where 
f(x ) is a monotone function. As we shall see, the same technique can be used for 
lattice points. Then we shall turn to inequalities. The most useful inequality in 
analysis is the Cauchy-Schwarz inequality. 

Following that, we shall prove some simple estimates which we shall use re¬ 
peatedly. Finally, we shall give a very simple evaluation of the Gauss integral. 
From the Gauss integral we shall derive the Hecke integral. 

2. Comparison of the Sum f(k) with the Integral fo f(x)dx 

In this section, we want to obtain some useful elementary estimates. We will 
use such estimates throughout the book. 

Consider the following figure. 
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A simple comparison of areas yields the estimates 


a) 


N N-l 

£/(*)< f" f(x)dx f(k). 

k=l k =0 


From these inequalities, it is easy to obtain asymptotic results, as will be indi¬ 
cated in the exercises. 

Exercises 

1. Derive a similar result for the case where f(k) is first increasing and then 
decreasing. 

2. Derive the two-dimensional version of these inequalities. 

3. Consider the case where the function of two variables has the form f(xy) 
and hence derive estimates for the sum 2^ =1 d(ri)f(n) where d(n) is the number 
of divisors of n. 

4. Consider the case where the function of two variables has the form f(x 2 + y 2 ) 
and thus obtain estimates for 2^ =1 r(n)f(n) where r(ri) is the number of repre¬ 
sentations of n as the sum of two squares. 

3. Lattice Points 

The same techniques can be used to count the number of lattice points in a 
region. This is part of the geometry of numbers which we shall not enter into 
here. 

Exercises 

1. Show that the number of lattice points in the region 1 < xy < n is 
[nfx] where [x] is the greatest integer less than or equal tox. 

2. Show that the number of lattice points in the quadrant of the circle x 2 

+ y 2 = n is given by [ n 2 - k 2 ] . 

3. For a convex curve, compare the number of lattice points with the area. 

4. Consider the triangular region determined by the lines x = Q,y = 0, ax + by 
= 1. Count the lattice points in this region in two ways and equate. 

5. Consider the tetrahedral region determined by the planes x = 0, y = 0, 
z = 0, ax + by + cz = 1. Count the lattice points in this region in three ways 
and equate. 



1.4 Estimates 


3 


6. Do the same for a quadrant of an ellipse. 

7. What meanings can be ascribed to the statement “The probability that two 
integers are relatively prime is 6/7T 2 ?.” (Lehmer, Bellman-Shapiro). 

4. Convergence 

We can do much better in many cases, let us write 


(i) 


N 

X PN +1 

f(k) -J f(x)dx 


k= 0 ' 


Let us now use the elementary result 
(2) m-m=Jj'(y)dy. 

We thus have 


N v N k 

(3) IrXi^T a f'(y)dy. 

k=l \ ~ * k=l k~l 

In many important cases, the function f'(y) is monotone. Consequently, we 
see that if $of*(y)dy converges, we have the limit. 


(4) 



Exercises 

1. Show that e* - 1 < xe*, if x > 0. 

2. Show that 1 - e~ x >0, if jc > 0. 

3. Show that 1 +x+ • * + =(1 -x« +1 )/0 -x). 

4. Hence, show that 
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5. Hence, by integrating, show that 

log —— = x + —+••• + f X —— dy. 
1-x 2 Jo 1 -y 

6. Using the representation 


iog(i +*>= r 

Jo 1 + y 


show that (1 + xjn) n is monotonically increasing in n. 

7. Assume that f(h) is monotone decreasing and convex. Obtain a better upper 
bound using the trapezoid formed by connecting k and k + 1. 

8. Obtain a lower bound by using the trapezoid formed by talcing the tangent 
at {k, f(k )). 

9. Obtain a lower bound by using the trapezoid formed by taking the tangent 
at (k + 1/2,/(£ + 1/2)). 

5. Euler’s Constant 


The most important example is Euler’s constant. We have 

(1) lim ( V - - logTV) = y. 

\ n ! 

The values of y is 0.5771 * • *. Most likely, y is transcendental, but it has not 
even been shown that y is irrational. The question of the arithmetic nature of 
y is one of the outstanding problems of the theory of numbers. 

We shall encounter the Euler constant in Chap. 4. 


Exercises 


1. Find the 

2. Obtain a 

3. Find the 


asymptotic behavior of n a e~ nx as x 0. 
two-dimensional generalization of the formula above, 
asymptotic behavior of 
' N 


I 


1 

nlogn 
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6. Inequalities 

In this section, we will derive the Cauchy-Schwarz inequality. 

We begin with the fundamental result 

(1) {a-b? >0. 

Multiplying out, we obtain 

(2) a 1 + b 2 > lab . 

Inwards, the arithmetic mean of two positive quantities is greater than the geo¬ 
metric mean. 

We now take the fundamental inequality in (2) and set 


(3) 



Adding, we obtain the famous inequality of Cauchy-Schwarz 


(4) 



This is one of the most useful inequalities of analysis. 


Exercises 

1. Using (2) above obtain the inequality a\ + a\ + a% + a% > 4a l a 2 a 3 a 4 • 

2. Continuing in this fashion, obtain the inequality for any power of two. 

3. By specialization obtain the famous arithmetic geometric mean inequality: 
The arithmetic mean of n nonnegative quantities is always greater than or equal 
to the geometric mean. 
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4. Show that the Cauchy-Schwarz inequality may be written 


N 


n =1 


i 1/2 


TV 


max 


TV 


2 * 




where the max 
is taken over; 


2 «-> 
n =1 


5. From the foregoing representation derive the Minkowski inequality 


( TV ' \ 1/2 / TV \ 1/2 /TV \ 1/2 

5 + (i«j 

6. Prove the arithmetic mean-geometric mean inequality for the case n - 4 
by showing that the difference is nonnegative. 

7. Establish the corresponding result for any power of 2. 

8. By specialization, obtain a corresponding result for general n. 


7. Continuous Version 

The continuous version of the Cauchy-Schwarz inequality is also very useful. 
We begin as before with the inequality 

(1) (f-g) 2 > 0. 

Multiplying out, and integrating we have 



This shows that the product fg is integrable whenever/ 2 and# 2 are. 
Using 



( 3 ) 
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we obtain 

{pf(/*Y 

8. A Technique of Shnirelman 

From the Cauchy-Schwarz inequality we have 



( 1 ) 



Solving, we have 



By a judicious choice of the a n , we can obtain a lower bound for the number 
of elements of a set S. 

9. Summation by Parts 

A most powerful technique for the estimation of integrals is integration by 
parts. A similar technique exists for sums, summation by parts. 

(1) ?n s n ~ tn-\ s n-\ ~ ^ni s n ~ s n -1) + s n -1 (fn ~^n- 1)* 

Summing over n, we have the formula 

N N 

(2) S N^N “ s 0^0 = tn(. s n “ s «-l) 

n ~1 n -1 




Exercises 

1. Show that we can estimate ^ n a n if we have the order of magnitude 


N 
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2. When can we go in the other direction? 

3. Obtain the order of magnitude of d(n). 

4 . Obtain the order of magnitude of 2^ =1 r(n). 

5. Use integration by parts f^ n \y)(x - y) n dy and thus obtain a representa¬ 
tion for a Taylor series with a remainder. 

6. Obtain a remainder estimate for the function e x . 

7. Obtain an estimate for the remainder of 1/(1 -x) a . 

8. Show that there is a similar formula for Appell polynomials defined by the 
relation p n (x) = p n - X {x). 

9. From the representation 

rx+y 

(x +y) s -x s = sj z s ~ l dz. 


derive the estimate 


--—-< yx° 1 , where o is the real part of s. 

s 

for 

0 <x,y, o. 

10 . Show that integration by parts may be obtained from summation by parts 
by going to a suitable limit. 


10. The Gaussian Integral 

One of the most important integrals in analysis is the Gaussian integral. There 
are many ways of evaluating it, and we will meet it later on in Chap. 4. 

We have 

(1) r r e <x2 * y2) dxdy = r r e-'rdrde = rr/4, 

J 0 d 0 J 0 Jo 

upon using polar coordinates. 

Thus, we have the famous formula 
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Exercises 

1. Show that the change of variable x =x' + ia is legitimate by considering the 
contour 


■R + ia 



R + ia 





-R 

R 


and letting R approach infinity. 


2. In this way evaluate the integrals poo e~ x cos xydx, poo e~ x s inxydx. 

11. The Hecke Integral 
In the Gaussian integral make the change of variable 

(1) x=y-w/y. 

This has the effect of mapping the interval (-°°, °°) into the interval (0, °°) 
in a one-to-one fashion. 

We see that we have 

(2) dx=dy(l+w/y 2 ). 

Making the change of variable y = 1//, we see that the second term yields 
essentially the same integral. We can obtain the same integral by making the 
change of variable y = wy . 

Combining the foregoing results, we have the expression 


( 3 ) 



e -y-wfy 



= e' 2 ^. 


This integral was used by Hecke to obtain many identities in the analytic 
theory of numbers. 
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Exercises 

1. Show that the above integral holds if we merely assume that the real part of 
w is positive. 

2. Prove that the function defined by the Hecke integral satisfies a linear 
second-order differential equation. Hence, evaluate the integral. 
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Miscellaneous Exercises 

1. Establish the identity of Lagrange 

]£«? = 1/2 £ (a,*,-a,h,) 2 . 

i i ' i / (/ 

2. What is an integral analog of this identity? 

3. Show that this identity may be put in determinantal form and thus obtain 
a higher-dimensional generalization. 

4. Show that if u n+m < u„ + u m , u„ > 0, then u n/m approaches a limit 
(Fekete). 

5. Can the condition u n > 0 be relaxed? 

6. Can the preceding lemma regeneralize to matrices where the relation A 
< B means that B-A is nonnegative definite? 
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2. Transforms 


1. Introduction 

In this chapter, we want to study some transform methods. These methods 
are very important in analytic number theory and throughout analysis in general. 

In Sec. 2, we study Fourier series. In Sec. 3, we will give the fundamental 
theorem of Fejer. In Sec. 4, we use this result to establish the result of Parseval. 
In Sec. 5, we derive a very important Fourier series, which occurs in many parts 
of analytic number theory. 

As a natural generalization of Fourier series we have the Fourier integral. The 
most important result here is the Parseval-Plancherel theorem. 

We turn next to the Laplace transform. In this and the succeeding section we 
study the inversion formula and the formula for the product. We turn next to 
the Mellin transform. Again, we are interested in the inversion formula and the 
formula for the product. Using the Laplace transform, we derive the fundamental 
identity of Lipshitz, which plays an important role in analytic number theory. 

The theory of congruences naturally leads to the study of finite Fourier trans¬ 
forms. In Chap. 12, we shall show that the Laplace transform plays a funda¬ 
mental role in Tauberian theory. It may also be used to evaluate various integrals 
and to derive various identities, as we show in the exercises. The Mellin trans¬ 
form plays a fundamental role in the study of Dirichlet series. 

2. Fourier Series 

The technique of Fourier series is one of the most powerful in analysis. Here, 
we shall only sketch the outlines of the method. 
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We begin with the fundamental orthogonality relation, 


0) 



Anx -imx _ 


dx = 0, 


n=£m, 


= 2i r, m = n. 


Using this relation, we have 


(2) /(*) = ^ 

n=~°° 

we have 


(3) a n = J_ /* flx)e mx dx, 

2tt J - 7 r 

provided that the series is absolutely convergent. 

We can define the Fourier coefficients by means of the relation in (3). 
We write 

n~oo 

(4) fix)- £ 

n——°° 


The question arises as to when the series on the right actually converges. 
Let us write 


s N 


n=N 

= J a * einx 


n=-N 

n=N 


- 1 (u: 

n = -N ' 




-ihf>i%^ e * y ) dy 

i r* 

= — / /OX- )dy- 

2lt J -n 


( 5 ) 
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The convergence question for Fourier series is delicate because the kernel is 
not absolutely integrable. Consequently, some condition must be imposed upon 
f(y). The simplest condition, and the one frequently used in practice, is due to 
Dirichlet. If f(y) is a bounded variation, then we have convergence. In most 
applications, that means that f(y ) is piecewise differentiable. Then we have con¬ 
vergence at the points of differentiability. 

Exercises 


1. If {<p„(jc)}is complete and orthonormal, by which we mean 



'Pn(x)'fm(x)dx = 0 , 


m^n. 


we have a correspondence 


= 1, m = n, 


f~ 

n~ 1 

where 

a„ = J n f{x)y n (x)dx. 

The coefficients are well defined by virtue of the Cauchy-Schwarz inequality if 
feL\ 

2. If the sequence is almost-orthogonal, by which we mean 

/ " 'Pm(x)'4,(x)dx = a mn 

-7T 


where 


^ a mn x m x n <k^X 2 n 
m,n=Q n- 0 

for some constant k, we also have a correspondence 


/' 




n=1 


[See 
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Bellman, R., Almost-Orthogonal Series , Bull. Amer. Math. Soc. 50 (1944), 

517-519. 

Almost orthogonal series are used in sieve methods.] 

3. Show that 


fe« 


2 a " e 


for any integer a. 

4. Hence, we have a correspondence 

The infinite matrix is called a Toeplitz matrix. 


5. 


lim f sinnx ( —— — = 0. 

n ^°° J -rr \sinx x) 


6. Hence, evaluate 



sinx 

x 


dx. 


3. Fejer’s Theorem 

In this section we will establish the famous theorem of Fejer. 

This is consistent with the general theme of the book. It may be quite diffi¬ 
cult to obtain the behavior at an individual value. However, the behavior of a 
suitable average, or mean, will be quite smooth. 

Let us define 

(1) So+S. +•••+ _% 

N+l 

the arithmetic mean of the partial sums. 

It follows that 

(2) On = r f(.y)(- ■ 0 dy. 

J -7r 

The important thing is that the kernel is now nonnegative. This nonnegativity 
makes the discussion of o n very easy. We thus have the famous result of Fejer: 
If/OO is continuous, o n converges uniformly throughout the entire interval. 
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Exercises 


1. Show that 


o N 


k=N 

■I 1 , 


k=-N 


N+ 1 


|fc| a k e ikx . 


2. Show that every continuous function over an interval may be uniformly 
approximated by a trigonometric polynomial. 

3. Show that every continuous function over an interval may be uniformly 
approximated by a polynomial. 

4. Show that every continuous function over an interval may be uniformly 
approximated by a function with an absolutely convergent Fourier series. 

5. Show that if /is continuous the condition 


/: 


f(x)e mx dx = 0, n = 0, ± 1, 


implies that fix) is identically 0. 

6. Show that the condition 


/: 


f(x)x n dx = 0, 


« = 0 , 


implies that f(x ) is identically 0 if/is integrable. 

7. Use Fejer’s kernel to establish the value of 

8. If the series D£ =1 a n has its arithmetic mean of partial sums convergent, and 
if 2£ =1 ka k = o(n), then the original series is convergent. (This is the original 
result of Tauber which stimulated Hardy and Littlewood in their research.) 

9. Show that if a periodic function is analytic, it has an absolutely convergent 
Fourier series, and all its derivatives may be calculated from the series. 


4. Parseval’s Theorem 


We have 

f" sldx= ^ a\. 


(1) 
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We wish to extend this result. For this purpose we shall use the foregoing re 
suit of Fejer. However, first, we have Bessel’s inequality. 


( 2 ) 



This asserts that the left side converges whenever the right side exists. To estab¬ 
lish this we begin with 


N 

(3) ^ a\ =J^f(x)s„(x)dx 

n=-N % ^ 

and use the Cauchy-Schwarz inequality. 

We have 

« =/>>*■ 

k=-N \ 

We know that as N -+ 00 , on the left, because of the convergence of the series 
we see that the limit is « a*. On the right, because of the uniform con¬ 
vergence we have J\ f 2 (x)dx. 

We thus have the formula of Parseval 

(5) S al= h5l , f2( ~ x)dx - 


vahd for any continuous function f(x). 

Exercise 

1. Obtain Bessel’s inequality by determining the minimum of 

f (f 2 - Cl 01 - • • • -c„<p„) 2 dt. 

J o 

5. Fourier Series for the Fractional Part 


Consider the function x - [x]. This represents the fractional part of x, and is 
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an important number-theoretic function. We see that it is periodic of period one. 
Hence, we have the Fourier series, 


(i) 


oo 



«=l 


as an easy calculation shows. 

In the following, this series will be used for several purposes. 

Exercises 

I. Use Parseval’s theorem to evaluate 



m n 

n—\ 


2. In the same fashion evaluate E^ =1 1 /n 2k . 

6. Finite Fourier Series 

In many parts of number theory, we meet a sequence which is periodic. It is 
convenient to treat such sequences by a finite Fourier series. We write 

N 

( 1 ) /(«)= 2 a k e 2 ™ k/N . 

k=-N 

There is now no question of convergence, nor of the existence of Parseval’s 
relation. We shall not consider finite Fourier series because they are intimately 
connected with trigonometric sums, which we have agreed not to treat. 

Exercises 

1. Iff(x) has an absolutely convergence Fourier series we have 

oo 

/(*/ N ) = 2 a r e 2 ™ k l N 

r =-co 

oo 

= X' 2n,>fc/yv (^ +a ^ + -)- 

r--N 
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2. Find the finite Fourier series for f(k/]sf)e 2Ulka ^ N . 

7. The Fourier Transform 

Let us now consider the case when the interval is infinite. If we proceed 
formally, we see that the series becomes an integral. This procedure may be 
made rigorous. However, it is easier to begin with the expression for the infinite 
integral 

(1) F{y) = f f(x)e‘ xy dx. 

J —oo 

This expression is called the Fourier integral. We are interested here in two 
features, the inversion formula and the Parseval relation. In order to make this 
symmetrical, the integral is often written 

(2) F(y) = -L f f(x)e lxy dx. 

\J2lT J 

Then we have the elegant inversion formula 

(3) A*) = -A= f F(y)e~ ixy dy. 

V 27T j 

8. Parseval-Plancherel Relation 

We have 

(1) r\F{y)\ 2 dy = \f(x)\ 2 dx. 

J -oo J -OO 


We will employ this result in Chap. 12. A proof of this fundamental result 
may be obtained in various ways as the book by Titchmarsh shows. 

9. The Laplace Transform 

We write 


L(f) = F(s) = [~ e~ st f(t)dt, t > 0 . 

Jo 


(1) 
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We may regard this as a continuous version of a generating function. We have 
the inversion formula 

( 2 ) f e st F{s)ds. 

2m Jc 

Here C is a contour which lies to the right of the singularities of F. We shall 
take C to be a straight line. This inversion formula allows us to use the theory of 
a complex variable. 

We may regard the Laplace transform as a special case of the Fourier trans¬ 
form, or conversely. 

Exercises 

1. Use the Parseval-Plancherel formula to establish the result 

r^>*. 

Jo Jo S + IX 

where now F is the Fourier transform of f. 

2. Use the result (1 + ix/n) n -+ e lx to obtain the post-Widder inversion formula. 
(Other proofs of this inversion formula will be found in the books by Doetsch 
and Widder cited in the bibliography.) 

3. Consider the maximum transform defined by the relation max^o f(x)e~ xy 
= F(y). Show under appropriate conditions that we have /( x) = voin y> QF{y)e xy . 

See, 

Bellman, R. and W. Karush, On a New Functional Transform in Analysis: The 
Maximum Transform, Bull. Amer. Math. Soc. 67 (1961), 501-503. 

-, Mathematical Programming and the Maximum Transform, J. Soc. Indust. 

Appl. Math. 10 (1962), 550-567. 

-, On the Maximum Transform and Semigroups of Transformations, Bull. 

Amer. Math. Soc. 68 (1962), 516-518. 

-, On the Maximum Transform, J. Math. Anal. Appl. 6 (1963), 67-74. 

-, Functional Equations in the Theory of Dynamic Programming-Xll : An 

Application of the Maximum Transform, J. Math. Anal. Appl. 6 (1963), 
pp. 155-157. 

4. Consider the function defined by the relation 
h(x)= max (f(x 1 )+g(x 2 )). 

X j ■ X 2 X 

ma xh(x) = max/(x) + max^(x). 

x>0 x>0 K x>0 


Show that 
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(This relation plays an important role in some problems of dynamic program¬ 
ming. See 

Bellman, R. and S. Dreyfus, Applied Dynamic Programming, Princeton Uni¬ 
versity Press, Princeton, NJ, 1962. 

10. The Product Formula 

Consider the function 

(1) h(x) = J*f (x - y)g(y)dy. 


Then, we have the fundamental result of Borel 

(2) L(h) = L(f)L (g). 

The function h(x ) is usually called the convolution of /and#. The validity of 
Borel’s formula may be established under various conditions on/and#. 

11. The Identity of Lipshitz 

Let us now derive a fundamental identity. This identity plays an important 
role in the theory of numbers, as we will show. We begin with the Fourier expan¬ 
sion of x - [x] derived above. 

We have 

a, 

x n~ 1 ' 


On the left, we write 


( 2 ) 



e~ sx (x - ri)dx; 


performing this integration, and interchanging the summation and integration on 
the right, the legitimacy we will leave to the reader, we obtain the identity 


1 

1 -e* 


.i + i 


1 

n =1 


2 s 

s 2 + 4n 2 n 2 


0. 


(3) 
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Exercises 

1. Obtain the partial fraction expansion of 1/(1 - x rt ). (An interesting gen¬ 
eralization of partial fractions was obtained by Jacobi. See Jacobi, C. G., Ges- 
amnaelte Werke. 

2. Using this result, with a suitable choice of x and letting n approach infin¬ 
ity, derive the foregoing identity. 

12. The Mellin Transform 

Let us now introduce the Mellin transform. We write 

(1) M(f) = F(s) = [~ 

J 0 

We have the inversion formula 

( 2 ) m = f f°F{s)ds, 

2m J c 

where again C is a contour to the right of the singularities F(s). We shall employ 
this formula extensively. 

13. The Product Formula 

Consider the function 

a) h ( X )= r mg(x/t)dt . 

Jot 


Then we have formally, 

(2) M(h)=M(f)M(gl 


This result may be established under various conditions on /and#. It enables 
us to evaluate many integrals as we shall see. 
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Miscellaneous Exercises 

1. Show that if we consider a double Fourier series and if the function is of 
the form f(x\ + x\ ), then the Fourier transform is of the same form. 

2. Obtain in this way the Bessel transform. 

2 

3. Use the Laplace transform on the integral e~ tx dx and thus evaluate the 
Gaussian integral. 

4. Show that the Fourier expansion may be obtained from integral equation 
theory by considering the equation 

y *2ir 

f(x -y)<t>(y)dy. 
o 
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This is one of the fundamental identities of analytic number theory. It has been 
generalized by Hecke and Siegel. 

This identity was first pointed out by Lipshitz. As we shall indicate in Chap. 6, 
there are important generalizations. 
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This result can be obtained from general results in the theory of the complex 
variable, the expansion of Mittag-Leffler. 

Section 12 

The Mellin transform can be discussed as a special case of either the Laplace 
transform or the Fourier integral. See the book above by Titchmarsh for a dis¬ 
cussion of rigorous aspects of this integral. 




3. Congruences 


1. Introduction 

In this Chapter we want to introduce the residue notation following Gauss. 
We shall recall various facts which we shall use in the following. This is not 
intended to be an introduction to congruences, rather a review. 

In the second section, we introduce the congruence notation. What makes this 
notation so powerful is the fact that the congruences can be manipulated like 
ordinary equations. In Sec. 3, we prove the fundamental theorem of Fermat. 
With the aid of this relation in Sec. 4 we can treat the linear equation, and give a 
simple representation for the solution. In Sec. 5, we discuss some facts about the 
quadratic equations. Here we say a word or two about the famous quadratic 
reciprocity law. In Sec. 6, we turn to the general polynomial congruence. Here 
algebraic number theory is required, and we shall just state the results we need. 

Finally, we point out the connection between congruences and trigonometric 
sums. 

2. The Congruence Notation 

We write 

(1) a = b(p), 

read “a is congruent to b modulo p”. 

This means that a - b is divisible by p. Congruences may be manipulated like 
ordinary equations. 
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Exercises 

1. Show that a necessary and sufficient condition that the congruence 

'YfijXj = 0(p) 

/ 

is that the determinant = 0(p). 

2. If this condition is satisfied, what is the solution? 

3. Show that the polynomial determined by the conditions 

p-i 

^ k r p(k) = 0 (p), 0 < r < p - 1, 

o 


is uniquely determined. 

4. Show that these polynomials are orthogonal. 

5. Show that these polynomials satisfy a 3-term recurrence relation and find 
this recurrence relation. (These polynomials are analogs of the Legendre poly¬ 
nomials.) 

6. Consider the nonzero residues modulo p. Show that by means of the Eu¬ 
clidean algorithm we can determine a reciprocal of each nonzero residue. 

7. Show that the nonzero residues form a group, if p is a prime. 

3. Fermat’s Theorem 

In this section we want to establish the fundamental result of Fermat. Natu¬ 
rally, there are many proofs of this result, and we shall indicate another one in 
the exercises. 

Let us begin with a relation 

(1) afi = bj(p). 

Here a if a, and bj are residues modulo p. We take initially p to be a prime. 
Here a i9 a , and bj are nonzero residues. 

If the di are different the bj must be different. This follows from the fact that 
there can be no zero divisors modulo p. We now take the relation in (1) and 
multiply over all the residues. The product on the left must be the same as the 
product on the right. The result is 


( 2 ) 


a p ' x = lO), 
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the famous theorem of Fermat. 

We can go a bit further. The relation in (2) is equivalent to saying that 

( 3 ) cf' x =\+kp. 

Raising both sides to the pth power, we have 

(4) a*>-i> = i(p*). 

Proceeding inductably, we have 

(5) aP*<P-i) = l(/,* +1 ). 

It follows from the unique factorization theorem that we have the extended 
form of the theorem of Fermat 

(6) a Hn) = l(n). 

Here, <j>{n) is the Euler 0 function studied in Chap. 8. 

Exercises 

1. Consider the set of powers 1, a, a 2 , • • ■. Show that two powers must be the 
same and thus we have the congruence 

c r* = l(p), where k may depend on a. 

2. Show that if a? = 1 (p), then d\(p - 1). 

3. Show that the nonzero residue modulo p form a group. 

4. Derive Exercise 1 from this fact. 

5. Consider the two by two matrices whose elements are the residues modulo 
p. Show that the nonsingular matrices form a group. 

6. Show that the matrices whose determinant is one form a group. 

7. What are the orders of these two groups? 

8. Obtain a Fermat theorem for matrices. 

9. Consider the recurrence relation 

u n+l = a ll u n +a 12 v n(P)> 
v n+ 1 =a 2 iu„ +a 22 v n (p). 
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Show that we have a Fermat theorem. 

10. Consider the recurrence relation 

u n+2 =au n +i +bu n (p) . 

Show that the sequence (u n ) is periodic and satisfies a Fermat relation. 

11. Prove the result of Lagrange that the order of a subgroup must be a divisor 
of the order of the group. 

12. Consider the set 1, a, a 1 , • • •, a p ~ 2 . Since there are p - 1 numbers, there is 
a possibility that all the residues are given this way for a suitable a. If such an a 
exists, it is called a primitive root. How many primitive roots are there? 

13. Show that if 

/(z) = (az + b)l(cz + d), 
then the iteration of /(z) has the same form. 

14. Consider the case where a, b , c, and z are residues modulo p . Show that the 
iteration is equivalent to matrix multiplication. 

15. Hence, derive a Fermat’s theorem for the iterates of/(z). 

4. The Linear Equation 

Consider the equation 

(1) an+b=0(p). 

By a simple change of variable, we may take the equation in the form 

(2) an = b(p). 

Since there are no zero divisors, there is at most one solution. Using Fermat’s 
theorem, we may exhibit a solution in the form 

(3) n=a p ~ 2 b(p). 

Here, as usual, p denotes a prime. 

Exercises 

1. What is a number of solutions of an = b(p) for 1 < n <x? 

2. Consider the equation AX = B(p ), where A, X, and B are matrices. 




3.6 Congruences 


35 


5. The Quadratic Equation 

Consider the equation 

(1) an 2 +bn +c = 0(p). 

By the usual change of variable, we may take the equation in the simpler form 

(2) n 2 =d(p). 

The first fact we observe is that this equation is not always soluble. Thus, for 
example, the equation 

(3) 2(5), 

is not soluble. Thus, we have the interesting question of what equations are solu¬ 
ble. If the equation in (2) is soluble, we call d a quadratic residue. If the equa¬ 
tion is not soluble we call d a quadratic nonresidue. 

It turns out that there is a relation between the solubility of (2) and 

(4) « 2 =p(d). 

This is the substance of the famous law of quadratic reciprocity, conjectured 
by Euler on the basis of examining a large number of cases, and proved by Gauss. 
In Chap. 6, we give an extensive generalization of this law involving Gauss sums. 

Exercises 

1. Prove that the quadratic residues form a group. 

2. If the congruence n 2 = a(p ) has a solution then ( p - 1 )/a 2 = 1(/?). 

3. Is the converse true? 

4. Is the product of two quadratic nonresidues always a non-residue? 

6. The General Polynomial Equation 

Let us now turn to the general polynomial equation 

( 1 ) m = 0(p). 

This equation needs algebraic number theory for its study. Let us cite the 
following results which we will use. Let p(p) denote the number of solutions of 
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this congruence. Then we have 


( 2 ) 



n~\ 


log N 


as N becomes infinite. We also have 

(3) V ~ log log N. 

P 

p<N 

Here we assume that/(x) is irreducible. 

Here, the summation is over the primes less than or equal to N. For the qua¬ 
dratic case, these relations may be established using the Legendre symbol. 

Exercises 

1. Show that every polynomial may be considered the characteristic polyno¬ 
mial of a matrix. 

2. Show that p(nm) = p(n)p(m), if (n, m)= 1. 

3. Show that we have the factorization 

nP - 1 s(w- 1) (n -2) •(« -p + 1 )(p). 


7. Connection with Trigonometric Sums 

We begin with the relation 

§ If pXr, 

k-0 P 

= P if P\r. 

We thus have a method for converting a congruence into the evaluation of a 
trigonometric sum. Thus, for example, we have the expression 


( 2 ) 


22 


n=Q k-0 


c 2mkf(n) 


= PP(«). 


P 
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= 0(p). 

i 

5. Show that the congruence f(n) = 0 (p) is equivalent to a congruence of the 
above form, and thus obtain a representation for p(n). (This method was sketched 
in, 

R. Bellman, A Note on the Solution of Polynomial Congruences , Boll. Un. Mat. 
Ital. 19 (1964), 60-63.) 

6. Show that 

= 0(y/p logp). (Polya) 
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Miscellaneous Exercises 


1. Show that two polynomials of the form x 2 + y 2 have a product of the 
same form. 


2. Show that this relation is equivalent to the addition formulas for the cosine 
and sine. 


3. Show that we have the factorization x 3 + y 3 + z 3 - 3 xyz = (jc + y + x) 
(x 2 +y 2 +z 2 - xy - xz - yz). 

Even more dramatic, as H. N. Shapiro points out, 2(x 3 + y 3 + z 3 - 3xyz) = 
(x + y + z) ((x - y) 2 + (y - z) 2 + (z - x) 2 ), which shows that the real zeros of 
x 3 + y 3 + z 3 - 3 xyz are just the plane x +y +z = 0 and the line x =y = z. 

4. Show that this relation can be used to prove that the cubic congruence an 3 
+ bn 2 + cn+d = 0(p)is equivalent to a simple cubic and quadratic congruence. 

5. Show that the congruence an 4 + bn 3 + cn 2 +dn+e = 0 (p) can be solved in 
terms of a simple fourth degree, third degree and quadratic congruence. 

6. Show that a similar result does not hold for the congruence of the fifth 
degree. (For this result, a knowledge of Galois theory is required.) 

7. Show that the result of Exercise 1 may be obtained by observing that 

x 2 +y 2 = (x + iy) (x - iy ). 


8. Show that the result of Exercise 1 may be obtained by observing that 


x 2 +y 2 


xy 

-yx 


9. Show that we have the correspondence 


10. Show that the product of two expressions of the form 

x 3 +y 3 +z 3 +3 xyz 
is again an expression of the same form. 

11. Show that the product of two expressions of the form x 2 - 2 y 2 is again 
of the same form. 

12. Given two polynomials of the form x 2 + y 2 , there is always a third such 
that the first is the product of the second and third. 

13. Derive Exercise 5 from the addition formula for the cosine and sine. 

14. Derive Exercise 5 from the representation x 2 +y 2 = (x + iy) (x - iy). 
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15. Use the representation 

x\ + x 2 + x 2 + x\ =(jCi + ix 2 +/x 3 +kx 4 )(Xi -ix 2 -jx 3 - far 4 ), 

where /, /, and k are quarternions to prove that the product of two polynomials 
which are the sum of four squares is always the sum of four squares. 

16. Derive Exercise 5 from the fact that x 2 + y 2 is the determinant of the 
matrix 


a 

17. Consider the congruence x 2 + y 2 = l(p). Associate with a solution of this 
congruence the complex number x + iy. Prove that the product of two such 
numbers yields another solution. 

18. Prove that the quotient of two such numbers yields another solution. 

19. Prove that these numbers form a group. 

20. What is the order of this group? 

21. Establish similar results for the congruence 

x\ +xj + x\ + X 4 = 1 (p). 

22. Consider the congruence xf + x\ + x\ = 1 (p). Show that we can obtain 
solutions by adding the conditions x 4 = 0(p). 

23. Show that no similar result exists for a sum of three squares. (A famous 
result of Hurwitz asserts that we only have a product formula for two, four, or 
eight squares. 

24. Consider the congruence y 2 = ax 3 + bx 2 + cx + d(p). Show that the addi¬ 
tion formula for Weierstrassian elliptic functions can be used to generate solutions. 

25. Consider the congruence y 2 = ax A + bx 3 + cx 2 + dx + e(p ). Show that the 
addition formula for Jacobian elliptic functions can be used to generate solu¬ 
tions. (The fact that elliptic functions enter into the consideration of these con¬ 
gruences is not surprising. Hasse showed that the Riemann hypothesis in fields 
of characteristic p can be treated by elliptic functions in the cubic case. Weil 
considered the general case using Abelian functions. 

26. Show how to generate solutions of the congruence x 2 - 2 y 2 = l(p). 

27. Show how to generate solutions of the congruence * 3 + y 3 +z 3 - 3 xyz = 
1 (P). 

28. Show that x n - 1 and x m - 1 have only the factor x - 1 in common if n 
and m are relatively prime. 
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29. Prove that the numbers c for which the congruence x 2 + y 2 = c(p) has a 
solution form a group. 

30. Prove that there exist values of p for which the congruence x 2 + y 2 = 
0(p) has a nontrivial solution. 

31. Consider that the quantities of the form a + ib, where a and b are not both 
simultaneously zero, form a group for some value of p and not for others. For 
what values of p do they form a group? 

32. Does there exist a polynomial g(n ) such that we have the congruence 
g{an + bn 2 ) = ag(n) (p)? 

33. Do we have a Fermat’s theorem for the iterates of an + bn 2 ? 

34. Does a primitive root exist for quantities of the form a + ibl 

35. Consider the 2 X 2 matrices where the elements are residues modulo p. 
Show that multiplication is not necessarily commutative. 

36. Show that the 2 X 2 matrices cannot be generated as powers of a fixed 
matrix. In other words, there is no analog for matrices of a primative root. 

38. Consider the general polynomial congruence. Multiply by the quantities 
x k f 0 < k < p - 2. Regard the resulting system as a set of equations for the quan- 
tites x k . Obtain in this way a necessary determinantal criterion for a solution. 

39. Is this condition sufficient? 
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4. The T Function 


1. Introduction 

In this chapter, we want to study the Euler T-function. 

This is one of the important functions of analysis. However, we shall not study 
it for itself. We are primarily interested in its application to the Riemann zeta 
function. 

We shall begin by establishing the fundamental functional equation. Then we 
shall provide an asymptotic estimate using the method of Laplace. Following 
that we shall establish some formulas using a representation due to Weierstrass. 
We shall mention some important extensions due to Siegel and Ingham. 

We have merely sketched the principle result. The reader is urged to fill in the 
details. 

2. The Basic Functional Equation 

The T function of Euler is defined by the equation below for R e(s) > 0 


( 1 ) 





Richard Bellman, Analytic Number Theory: An Introduction ISBN 0-805 3-0452-5 

Copyright © 1980 by The Benjamin/Cummings Publishing Company, Inc., Advanced Book 
Program. All rights reserved. No part of this publication may be reproduced, stored in a 
retrieval system, or transmitted, in any form or by any means, electronic mechanical photo¬ 
copying, recording, or otherwise, without the prior permission of the publisher. 


45 





46 


Analytic Number Theory: An Introduction 4.2 


An easy integration by parts yields the fundamental functional equation 

( 2 ) r(s+ l) = sr<». 

By means of this equation, we can easily obtain the analytic continuation 
throughout the entire complex plane. Furthermore, we obtain the fact that the 
function defined this way has simple poles and we can obtain the residue at each 
pole. 

The same integration by parts that yielded the fundamental functional equa¬ 
tion shows that we have 

(3) l\n + 1) = /i! 

for any integer n. Thus, the T function may be considered as a continuous 
interpolation to the factorial. 


Exercises 

1. Show that the analytic continuation may also be obtained by writing the 
integral as a sum 



e~ x dx + 



e’ x dx. 


2. Show that the quotient of any two solutions of the functional equation is 
periodic. 

3. Show that 

r'(s + i) _r'(s) + i 
r(s +1) r(s) s * 

4. Show that the integral / 0 °° e k ^ x x a dx may be expressed in terms of the T 
function. 

5. Show that the integrals 

r°° cosx. r°° sin* , , 

/ - d Xy I - ax where 0 < a < 1 

J o x a Jo x a 


may be expressed in terms of the T function by using the following contour. 
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6. Show that 


f°°e~ xy x s ~ 1 dx = 02- . 

Jo ./ 

3. r(s) as a Real Function 

Let us now briefly consider T (s) as a real function. 

It is obvious by differentiation that we have 

(1 ) r"(s) = J**° (log jc) 2 * 5 ' 1 e~ x dx. 

From this, it follows that T(s) is a convex function for real s. 

In the next section, we shall obtain the asymptotic behavior. 

Exercises 

1. Determine where T(s) has its minimum value. 

2. Show that the only logarithmically convex solution of the functional equa¬ 
tion is a multiple of the T function. (Artin) 

4. Laplace’s Asymptotic Evaluation of an Integral 

In this section, we follow a method of Laplace to obtain the asymptotic evalu¬ 
ation of T(s). 
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Let us recall a result due to Stirling. 

(1) n ! ~ n n e~ n \j2irn . 

This formula is very important in probability theory. To obtain this result and 
its generalization to the T function we use a method due to Laplace. 

Let us write the integral as 


( 2 ) r(n+s) = f°°x 

Jo 

If we now make the change of variable x = ny, this has the effect of making 
the integral 

(3) n n+1 j°°x s - l e- n V°*x-x)dx. 


This change of variable has the advantage of localizing the contribution of the 
integral at x = 1. We observe that the function log* - x has its maximum value 
at x = 1. Consequently, we expect that as n becomes larger, the contribution of 
the integral will center about x = 1. It remains to make the argument rigorous. 
To do this we write the integral as the sum of three parts. 


(4) 


r<*> f\-e ri+e n°° 

Jo Jo J 1-e J 1+e 


The choice of e will be made later. 

In the interval (1 - e, 1 + e), we use a Taylor expansion for the function 
x - logx. It has the form 

(5) x - logx = 1 + c 2 (x - l) 2 + • • •. 

We know that the radius of convergence will be determined by the first singu¬ 
larity of the function. That occurs at x = 0. Consequently, if e is chosen small 
enough, we know that the series converges. However, we do not want the full 
series. Rather, we want the first two terms and an error term. We know from 
exercises in Chap. 1 that we have 

(6) x - logx = 1 + c 2 (x - l) 2 + Oflx - ll) 3 . 


Consequently, we have 
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(7) f°°x s ' 1 e~ n(x - loeJc) dx= f 1+e x s ~ 1 e~ nll+C2(x ~ 1)2+0(x - 1)3l +/t„. 

Jo Jl-e 

We leave it to the reader to show that R n is of small order as n gets large. 

It remains to choose e. We want to choose e so that R n is small and the contri¬ 
bution from the (x - l) 3 is small. We can do both. A suitable choice of e, for 
example, is l/\/n. We can write 

(8) e -«(i-*) 3 = i +0 (n(l-x) 3 ) 

in the interval (1 - e, 1 + e), where the 0-term is uniformly bounded by 2 n 
(1 -x) 3 . 

The trick is now to expand the interval (1 - e, 1 + e) to the full interval (0, °°). 
We thus find that the asymptotic evaluation depends upon the integral 

/ oo 2 

^ 6* dx 

which we have evaluated in Chap. 1. This explains the appearance of the con¬ 
stant \/7r. 

We thus obtain the result 

(10) r(n +s)~ 2^, 

an interesting generalization of Stirling’s formula. 

Exercises 

1. Show that the T function cannot have any zeros for any real value. 

2. Obtain the asymptotic evaluation by using the functional equation for the 
derivative of the T function and Euler’s constant. 

5. The Representation of Weierstrass 

Let us now give an important representation of Weierstrass for the T function. 
We will use this representation to obtain many other results. 

This representation can be obtained from general theory. However, we will use 
the asymptotic evaluation. The basic idea is quite simple. We shall iterate the 
functional equation, use the asymptotic evaluation, and then group the terms 
suitably. 

Let us now go through this procedure. We have 
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(H 1 _ s ^ s(s + l)-•-(s + ri) 

r (s) r(s + 1) r(s + n + 1) 

upon iteration. We can write this 

r($) r(5 + /i + i) 

Unfortunately, the infinite product 


(3) 



diverges. Consequently, we introduce the convergence factor e^ n and write the 
product in the form 


oo 



This infinite product converges. Using the Euler evaluation, 


n 

(5) = lo gn+7 + °(^) 

k-1 

and the asymptotic form of the T function we obtain the representation of 
Weierstrass 


( 6 ) 



Exercise 

1. Show that we need only the asymptotic behavior of V(ri)/r(n + s). 
6. Half the Sine Function 


Using the Weierstrass representation, we have 
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a) —-— 

mn-s) l \\ n 2 ) 

n-\ 

sin ns 

~ -s -. 

IT 

Since the infinite product converges absolutely and uniformly, we can group 
the terms conveniently. 

Using the functional equation we obtain finally the formula 

(2) r<»r(i -*)=—?-. 

sin 7 rs 

Exercise 

1. Show that the foregoing result yields an evaluation of r(%). 

7. The Duplication Formula 


Consider the function 



Using the functional equation, we see that this function is periodic of period 
Vi. Using the asymptotic evaluation, we see that this periodic function is ac¬ 
tually a constant. We leave it for the reader to determine the constant. 

There are many ways of establishing this formula. Below we shall establish it 
in the manner of Euler using an infinite integral. It may also be established by 
using the results of a theory of a complex variable. 

We also leave it to the reader to use the Weierstrass representation to obtain 
this formula. 

Exercises 

1. Show that 


r(«s) 


is a simple function, and evaluate this function. (Gauss) 
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2. Obtain the duplication formula by using the Hecke integral and the Mellin 
transform. 

3. Use the result of Gauss and the Mellin transform to evaluate the integral 



for suitable values of a and b. 

8. The Method of Stationary Phase 

In Sec. 4, we considered the asymptotic behavior as s approaches infinity along 
the real axis. In this section, we want to consider the asymptotic behavior as 
s approaches infinity along the complex axis. 

The method again is simple. Using an appropriate contour, we show that the 
change of variable x = iy is legitimate. As s approaches infinity, which means as 
y approaches infinity, the rapid oscillation of the function makes the integral 
small. It is least small in the neighborhood of a stationary point. Consequently, 
we now repeat the procedure of Sec. 4. We will leave all the details to the 
reader. 

9. The Beta Function 


Let us write 

a) 



x m -\l -xf-'dx. 


This integral, the Beta function of Euler, is defined for 
(2) Re(m) > 0, Re(w)>0. 

It occurs in many parts of analysis. As we shall show, it can be easily evaluated 
in terms of the gamma function. 

We begin with the change of variable y = xj{ 1 - x). 

This transformation has the effect of mapping the interval (0, 1) into the 
(0, °°) in a 1-1 fashion. 

We then have to evaluate the integral 


(3) 



y m “1 

(1 +y) m+w 


dy. 


To evaluate this integral, we use the representation 
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(4) 


1 

(1 +y) m+n 



g-il+y^m+n-lfa 


Using this representation, and inverting the order of integration, we obtain the 
desired result 


(5) 

Exercises 


B(m,n) = 


r(m)rQi) 

T (m + n) 


1. Using the evaluation of the Beta function and the change of variable nx = 
y, obtain the asymptotic expression 

r(n)/r(n+s)~n- s . 

2. Evaluate the integral Iodx/\/x(l -x) by means of integration and the Beta 
function and thus evaluate TQA). 

3. Show that the integrals /J (sin 0) rt (cos S) n dd can be evaluated in terms of 
the Beta function. 

4. Show that the ffx m y n dxdy over the region x 2 + y 2 < 1 can be expressed 
in terms of the beta function. (Liouville) 

5. Find the area of the ellipse x 2 /a 2 +y 2 /b 2 = 1. 

6. Find the area of the circle and thus evaluate T(^). 

7. Find the volume of the ellipsoid x 2 ja 2 +y 2 /b 2 + z 2 jc 2 = 1. 


10. The Integral of Siegel 

An extensive generalization of the Euler integral exists, 

(1) f i^r 1 e +( A *>dA =(*••)• 

Ja> o 

Let us explain the notation. The integral is over positive definite matrices. 
The notation \A\ indicates the determinant of A. The expression Xx(AB)> is the 
sum of the elements along the main diagonal. 

The evaluation of this integral can be carried out inductively. In the bibliog- 
graphy and comments, we give references to some generalizations of this inte¬ 
gral, which forms an essential part of Siegel’s theory of matrix automorphic 
functions. An equivalent integral was found by Ingham in connection with 
mathematical statistics. 
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Miscellaneous Exercises 

1. Show that T(s) satisfies no polynomial equation. 

2. Show that r(s) satisfies no polynomial differential equation. (Holder) 

3. Prove that T(x) satisfies the recurrence relation 

r'(* + 2) = (2x - i)r\x +1) -*r'(*). 

4. By considering the quotient T*(x)/r(x + 1) obtain a formal continued 
fraction. 

5. In this way obtain a continued fraction for c. 

6. Does this continued fraction converge? (For questions of this type, see 
Perron, 0., Die Lehre von den Kettenbriichen , Chelsea, New York, 1950. 
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5. Riemann Zeta Function 


1. Introduction 

In this chapter, we will study some properties of the Riemann zeta function. 
This is a very interesting function and well merits study for itself. However, we 
are primarily interested in it as a generating function for the elementary arith¬ 
metic functions. 

In the bibliography at the end of the chapter, we will give references where 
further information can be obtained. 

2. Dirichlet Series 

The Riemann zeta function is defined by the series. 


(1) ?(s) = ^V s , Re(s) > 1 • 

n =1 

This series is the most famous example of a class of infinite series which are 
very important in analytic number theory, the Dirichlet series. The general form 
of a series of this type is 

oo 

(2) /(s) = ^a n n s . 

n =1 


It is clear that the sum of two series of this form is again a series of this form. 
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What is important is that the product of two series of this form is again a series 
of this form. We have 


(3) 

where 


f(s)g(s) = ^c n n s 

n=l 


( 4 ) c„ = ^ a k b, 

kl=n 

= ^ a k^n/k- 
kin 

Exercises 


1. Prove that 



2. Prove that the series on the right converges for Re(s) > 0. 


3. The Euler Product 


Let us now give the most important result concerning the zeta function. We 
begin with the observation that we have the unique factorization theorem. The 
analytic equivalent of this is the famous representation of Euler 


<■> iy- nKy 

n-\ p ' ' 

What is remarkable about this result is on the left we have summation over 
the integers, while on the right we have a product over the primes. This formula 
is the basis for all the work on the asymptotic number of primes. 

Exercises 


1. Obtain the Dirichlet series for log f(s). 

2. Obtain the Dirichlet series for f'(s)/f(s). 




5.5 Riemann Zeta Function 


59 


4. The Mobius Function 

In general, it is quite difficult to find a reciprocal of a Dirichlet series. If, how¬ 
ever, the Dirichlet series has an Euler product it is quite simple. Thus, we have 


« «■>"=§f=n(‘v> 

The function fi(n) is one of the most interesting functions of number theory. 
As we shall see, its behavior is intimately connected with the location of the 
complex zeros, of the zeta function. 

The function j u(ri) can be defined from this representation, or independently as 

m(«) = 1 , if n has an even number of simple prime factors, 

(2) = -1, if n has an odd number of simple prime factors, 

= 0, otherwise. 

If we use the obvious result 

(3) 

we have the fundamental relation 


( 4 ) ^ fx(/c) = 0 , 1 , 

Jc/n 

= 1, n = 1. 

Exercise 

1. Regarding as f'(s) times l/f(s), obtain an identity for the coeffi¬ 

cient of f'(s)/f(s). 

5. The Mobius Inversion Formula 

Using the properties of the Mobius function we obtain an interesting inver¬ 
sion formula. 

We have 

oo 

fix) = ^g(nx), 

n =1 
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(i) 

oo 

g(x) = £ n(n)f(nx). 

n=l 

This formula may be established under various conditions. Formerly, it follows 
from the properties of the Mobius function by rearrangement. 

6. The Squarefree Function 

If we square the mobius function, we obtain a very interesting counting func¬ 
tion. We have 

u 2 (n) = 1, if n is squarefree, 

( 1 ) 

= 0, otherwise. 

We shall study this function in more detail in Chap. 10. 

7. Multiplicative Functions 

It turns out that the elementary functions have a very important property. We 
introduce the concept of multiplicative functions. If 

(1) f(mn) = f(m)f(n) if (m, n) = 1, 

we say that the function /(n) is multiplicative. 

If 

(2) f(mn) = f(m) + /(«) if (m, n) = 1, 

we say that the function f(n) is additive. 

We have 

( 3 ) m=Hf(p a ) 

p 

if 


n 



P a . 


p 


(4) 
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It follows that we have a generalized Euler product. 


( 5 ) 



) 


Exercises 

1. Show that f(n k ) is multiplicative if f(ri) is. 

2. Show that if f(n) is additive and monotone, it must be a multiple of the 
logarithm. (Erdos) 

3. Show that if u(x) satisfies the equation u(x + y) = u(x) + u(y), and is dif¬ 
ferentiable, it is of the form u(x) = kx. (Cauchy) 

4. Show that the same conclusion holds if we merely assume that u(x ) is 
continuous. 

5. Show that the same conclusion holds if we merely assume that w(x) is 
integrable. 

6. Show that if we have the inequality | u(x + y) - u(x) - u(y)l < e, then a 
constant k exists such that \u(x) -kx\ <3e. 

7. Show that the constant 3 cannot be improved. 

8. Show that if u n > 0 and u m + n < u m + u n then ujn has a limit.(Fekete) 

9. If /„+! >f n - a n where 2^ =1 is convergent, then/ w converges. 

10. If/and ^ are multiplicative then the convolution of/and g is also. 

8. Analytic Continuation 

In this section we want to consider f(s) as an analytic function. 

The first task is to extend f(s) over the whole plane. This can be done in many 
ways. 

One way is to use the fundamental representation 


( 1 ) 


r 

Jo 


e~ x x s ~ x dx 
1 ~-e x ‘ 


We now write this integral as the sum of two parts. 
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The second term is an entire function of s. Let us then consider the first term. 
We have 




Consider the function 


( 4 ) 


e~ x 

1 - e~ x " x 


This function has its nearest singularity at 2m. Consequently, the power series 
converges for |x| < 2n. 

Integrating term by term, we obtain the expansion 


( 5 ) 


iW(s)= £ 

n—0 


S + n 5-1 


The last term comes from integration. This furnishes the desired analytic con¬ 
tinuation. The singularities at -n are cancelled by the zeros of 1 /r( 5 ). 


Exercises 


1. Show that 


2. Show that 



x-[x] 


dx. 



m 


for 0 < Re(5) < 1. 

3. Show that the formula 



4. Show that this formula gives the singularity at s = 1. Show that the other 
singularities do not exist. 
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9. The Functional Equation 

Let us now derive the famous functional equation for f(s). 

We shall employ the well-known functional equation for a theta function. We 
shall derive that equation in Chap. 6 and again in Chap. 7. Here we shall merely 
use it without proof. 

We use the representation 


0 ) 


' ' n~ 1 


-nt 


dt 


We write 

( 2 ) 


Next, we write 


(3) 


r* r* r 

JO JO J1 



In the interval (1, °°), we make the change of variable t = l/r, and use the 
transformation formula for the theta function. 



The final result is 

(5) = " (w) r {^y} f (1 - s). 

This shows that ifris/l) ?(s) is symmetrical about the line Vi + it. 
Exercises 


1. Let r(n ) denote the number of representations of n as a sum of two squares. 
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Establish a functional equation for 

Z, n ° ' 

n= 1 

2. Let r k (ri) denote the number of representations of n as the sum of k 2 . Es¬ 
tablish a functional equation for the associated Dirichlet series. 

10. The Riemann Hypothesis 

The famous Riemann hypothesis is that all the nontrivial zeros of f (s) are on 
the line Vi + it. 

Despite the best efforts of many mathematicians this has defied proof. As we 
shall see, it would have many interesting consequences for error terms. 

11. Asymptotic Behavior 

Fortunately, we have no need of a good point estimate of f (s) for our pur¬ 
poses. We can use a very crude estimate. To obtain this estimate, we use the 
series 



Our procedure is very simple. We break the sum into two parts, and estimate 
each part crudely. To obtain the required asymptotic estimate it remains to 
choose the way we decompose the sum. Let us write 


oo N ao 

® 2 - 2*1 

n =1 «=1 JV+1 


Choosing N to be [t ], we see we have the estimate 

(3) |f04 + Jf)l =0(^). 

We can obtain a better point estimate using a mean value estimate. Let us use 
the result 



dt = 0(TlogT). 


(4) 
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We shall obtain a more precise result in Chap. 12. As the exercises show, this 
result can be obtained in various ways. 

To obtain the point estimate, we use the fact that f (s) is an analytic function 
of s. Thus, we can employ the representation of Cauchy. 

( 5 ) m .± rm*L. 

2m Jr w - s 

It remains to choose the region R. We choose the following region 



V 2 + it 


Using this contour, we can estimate f (s) in terms of various mean values. To 
estimate the integral, we use the Cauchy-Schwarz inequality. We have 



The contribution from the other 3 sides of the rectangle are of lower order. 
The final result is 

(7) If (a + if) | = 0(log f 3/2 ) 

for a > x /i. 

Exercises 

1. Use the equation in (1) to establish the mean value result stated. 
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2. Estimate in a similar fashion the integral of f (s). 

3. Estimate in a similar fashion the derivative of f (s). 

4. What contours do we use if we want to estimate the function on the critical 
line. 

12. Uses of the f Function 

Once we have established the analytic properties of f (s) we can use it to obtain 
the behavior of many familiar functions of analysis. 

The basic idea is quite simple. We use the Mellin transform to obtain a repre¬ 
sentation, use the Mellin inversion formula, and then shift the contour. 

Let us give an example. We have 

( 1 ) r t s-i£tldt = r(sX(s + a). 

Jo n 

n =1 

From this follows 

(2) V1^! = i_ r r*r(s)f (s + a)ds. 

" n a 2m Jc 

n—\ 

If we shift the contour to the left, we obtain the representation 


(3) 



n=1 


_ r(i-«) 

t l-a 


f(a) + f f*r(s)?(s + d)ds. 

2m Jc ' 


Here C' is a contour which is a straight line to the left of s = 0. 


Exercises 

1. Find the behavior as t -+ 0 of 

2. Show that 

V I e — 

" n 1 - e~ nt 

n =1 


5) lo g(l - e ' nt ) = 2 


n=l 


-mnt 


m 


3. Hence, show that 
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where 


f[(l -e nt ) = e™ 
n=l 


/(0 = 



4. Hence, determine the asymptotic behavior as f ^ 0 of (1 - e~ nt ). 

5. Determine the asymptotic behavior as* -+ 1 of 2^ =1 jc"/( 1 -x”). 

6. Determine the asymptotic behavior as* and.y->0ofFQ „ =1 (1 -e~( mx+ny ^). 


13. Mean Values of the f Function 

We can also use the foregoing representation to obtain mean values. Remem¬ 
bering the dictum of Jacobi, we can use the foregoing representation of f (s). 
We have 


n=l 

where ^<a<l, a>0. 


14. The Prime Number Theorem 

Let us say a few words about the prime number theorem. 

Let us denote by 7t(x) the number of primes less than or equal to x. Analyti¬ 
cally, we have 


0) 


<*)= 2 i. 

p< X 


Tables showed that 
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Finally, in 1899, Hadamard and de la Vallee Poussin established this result 
rigorously. Their proof, which is classical, depended upon showing that the zeta 
function did not vanish on the line o = 1, and depended heavily upon the theory 
of a complex variable. 
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1. Write aXb to denote the convolution product. Show that 



aXb = bXa 


2. Show that aXbXc is unambiguously defined. 


3. Show that 



where 

HI 

tl 

t: 

-o 


4. Show that 

•t? 

ii 


where 

d («)=2 1 - 

k\n 


5. Prove that 

i 

/ \-i N 

n 4 

p<N \ n=l 


6. Hence, show that 

51 -^loglog'AT-d. 

" p 

p<N 


7. Show that 

logf(s) + c 2 > 51 -r 

p<iV ^ 


for any s greater than 1. 
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8. Choose 


s = 1 +- 


(logA0(log log AO 


and hence show that we have 


(1 +o(l))loglogjV> S' i-c 3 . 

" P 

p<N 


9. Hence show that 


— ~ log log N. 



10. Obtain conditions in f(n) so that we have a similar result for 

m 

p<N P 

1 1. We say the two integers m and n are relatively disjoint if they have no 
power of two in common in their representation. We shall use the notation 
[m, n] = 1. We say that a function is additively disjoint if it satisfies the func¬ 
tional equation 

f(m + n ) =/(m) +/(«, [m,n ]) = 1. 

Construct a theory of additively disjoint functions. 

12. Construct a theory of multiplicative functions for Gaussian integers. 

13. Is there an analog of the fundamental theorem of arithmetic for matrices? 
(Cauchy’s theorem was extended to matrices by Poincard, H., Sur les groupes 
continus, Trans. Cambridge Philos. Soc. 181(1899), 220-225. See the Bellman 
book Introduction to Matrix Analysis , McGraw-Hill, New York, 1960; 2nd 
ed., 1970.) 

14. Is there a theory of multiplicative functions for matrices? 

15. Is there an analog of the fundamental theorem of arithmetic for qua¬ 
ternions? (A theory of quaternion functions similar to the theory of functions 
of a complex variable was constructed by Fueter. See Fueter, R., Functions of 
a Hyper Complex Variable , University of Zurich, 1948-1949, reprinted by 
Argonne National Laboratory, 1959.) 

16. Is there a theory of multiplicative functions for quarternions? 
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17. Does the sene: 



n =1 


converge for t not equal to 0? 

18. Consider the functional equation 

Y(s + t)=Y(s)Y(t), Y(0) = l 

Show that if we assume that Y is differentiable, Y must have the form Y(t) 
- e At . Here Y and A are matrices. 

19. Show that the same conclusion holds if we assume that Y is continuous. 
(Polya) (See the book 

Bellman, R., Introduction to Matrix Analysis , McGraw-Hill, New York, 1960; 
2nd ed., 1970.) 

20. By using the relation f R $ (s)ds = 0, where R is a suitably chosen rectangle, 
derive a nontrivial estimate for /Jf (& + it) dt. 

21. We have the relation 

f «j)Ki-*)**o. 

Jr 

Can this relation, plus the functional equation be used to obtain the mean val¬ 
ues on the critical strip? 
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to Number Theory , Proc. Internat. Congress Math., 1954, Amsterdam, Vol. 
Ill, pp. 13-15, Noordhoff, Groningen; North-Holland, Amsterdam, 1956. 
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Soc. Transl. 19 (2) (1962), 47-85. 
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Galambos, J., Limit Distribution of Sums of (Dependent) Random Variables 
with Applications to Arithmetical Functions, A. Wahrscheinlichkeitstheorie 
und Verw. Gebiete 18 (1971), 261-270. 

Section 9 

A very interesting problem is whether this functional equation determines the 
function. As stated above, questions of this type were investigated in great de¬ 
tail by Hecke. See 

Bochner, S., On Riemann ’s Functional Equation with Multiple Gamma Factors, 
Ann. of Math. 67 (2) (1958), 29-41. 

Berndt, B. C., The Functional Equation of Some Dirichlet Series . II, Proc. 
Amer. Math. Soc. 31 (1972), 24-26. 

This was the second method used by Riemann to derive the functional equa¬ 
tion. See 

Riemann, B., Uber die Anzahl der Primzahlen unter einer Gegebenen Grosse, 
Werke 2, pp. 145-153, 1892. (Collected works of Bernhard Riemann, Dover, 
New York, 1953.) 

See 

Chandrasekharan, K. and H. Joris, Dirichlet Series with Functional Equations 
and Related Airthmetical Identities, Acta Arith. 24 (1973), 165-191. 

Section 10 

For a discussion of the Riemann hypothesis see the book by Titchmarsh cited 
above. With the aid of this hypothesis, many error terms can be improved. 

Hardy showed that there were infinitely many zeros on the critical line. For 
his method, see the book by Whittaker and Watson cited above. A. Selberg 
showed that a positive fraction of all zeros were on the critical line; for his 
method, see the book by Titchmarsh cited above. Levinson showed that more 
than one third of all the zeros were on the critical line using a different approach; 
see the paper, 

Levinson, N., “Complex Variables”. San-Fransisco, Holden-day, 1970. 
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Hilbert is reputed to have said, “If I woke after a thousand years, the first 
question I would ask is, ‘Has the Riemann hypothesis been settled yet?’” 
He also gave a talk in which he discussed three famous problems of mathe¬ 
matics: the Riemann hypothesis, Fermat’s last theorem, the ot? hypothesis. He 
said that so much had been done on the theory of entire functions that it was 
clear that the Riemann hypothesis would be settled in about twenty years at 
most. He said further that so much had been done in algebraic number theory 
that Fermat’s last theorem would also be settled soon. However, he added, the 
a* 3 hypothesis was forever beyond the power of mathematics. 

The hypothesis was done by Gelfand and Shneider independently in about 
fifteen years. 

See also, 

Good, I. J. and R. F. Churchhouse, The Riemann Hypothesis and Pseudoran¬ 
dom Features of theMobius Sequence, Math. Comp. 22 (1968), 857-861. 

The zeta function can be generalized in many ways. Perhaps most important 
is the zeta function in fields of characteristic p. This was done by Artin. By 
examination of particular cases, he was led to enunciate a Riemann hypothesis 
for these functions. This was proved in the cubic case by Hasse using elliptic 
functions. The general case was done by A. Weil using Abelian functions and 
algebraic geometry. 

The result has important applications to the estimation of trigonometric 
sums. 

Wintner, A., Random Factorizations and Riemann 's Hypothesis, Duke Math. J. 
11 (1944), 267-275. . 

Hejhal, D. A., A Remark of the Lindelof Hypothesis, Bull. Amer. Math. Soc. 
80(1974), 695-699. 

Section 11 

The pointwise asymptotic behavior of f(s) seems to be quite difficult. At 
present, the best results are obtained by the use of trigonometric sums, as dis¬ 
cussed in the paper, 

Vinogradov, I. M., The Method of Trigonometrical Sums in the Theory of Num¬ 
bers, Trav. Inst. Math. Stekloff 23 (1947), 109 (Russian). 

Section 14 

For an excellent account of the classical method for the proof of the prime 
number theorem, see the book 

Chandrakersan, K., The Elementary Arithmetic Functions, Springer, Berlin, 
1970. 



6. The Poisson Summation Formula 


1. Introduction 

In this chapter we shall discuss one of the most important tools of analytic 
number theory, the Poisson summation formula. We shall give three approaches 
to this fundamental result. As usual, the reason why we give different approaches 
is that each approach generalizes in a different direction, as we shall see. 

We shall discuss an approach using Fourier series. Then we shall present a 
method which depends upon Dirichlet series. Finally, we shall indicate how the 
Laplace transform may be used. 

We shall indicate at various places how this formula may be applied. 

2. Fourier Series 

Since our first method of proof of the Poisson summation formula will be based 
upon the theory of Fourier series, let us briefly digress to present the rudiments 
of this theory. Let us repeat a little of Chap. 2. 

Consider the set of exponential functions {e 2rrinx } 9 where n assumes the values 
n = 0, n = ±1, n -= ±2, and so on. These are clearly periodic functions of x of 
period one. It is plausible, from a number of physical considerations, to suspect 
that every continuous function of x which is periodic of period one can be 
represented as a linear combination of these particular functions. In other words, 
our physical surmise is that every periodic motion is a superposition of these 
simple periodic motions. 

Setting oo 

(1) /(*)= £ a n e 2 "™, 

n~-°° 
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how do we determine the coefficients? Fourier solved the problem by expanding 
both sides as power series in x and “solving” the resulting system of linear equa¬ 
tions for the a n . Following Euler, let us employ the orthogonality relation 

(2) [' e 2irimx ^ inx dx = 8f n . n . 

J 0 

Here 5 W _„ is the Kronecker delta symbol defined by the relations 5 m _* = 0, 
m =£ n f 5 m _ w = 1, m = n. It is reasonable then to suppose that if (1) holds, then 
the coefficient a n is determined by the simple formula 

(3) a n = f 1 f(x)e' 2nmx dx. 

Jo 

Oddly, despite the physical simplicity of the situation, its mathematical as¬ 
pects are both extraordinarily complex and subtle. The difficulty is essentially 
that the physical, intuitive concept of a continuous function is far too naive 
for the mathematical definition of a continuous function. Consequently, we cir¬ 
cumvent a number of irritating pitfalls by working backwards. We start with a 
continuous function /(x), (a type of function sufficiently general for our subse¬ 
quent purposes) and form the sequence of Fourier coefficients, {a„}, n - 0, 
±1, ±2, • • * , by means of the formula (3). Using these coefficients, we form 
the series 

oo 

(4) g(x)= X a n e^ inx 


where x is a real variable, lying in the open interval (0, 1). This series yields a 
new function £(x), defined where (4) converges. 

We then ask ourselves two questions: 

(a) Does (4) converge for allx, and, if not, for what values of x does it converge? 

(b) When (4) converges, is it equal to the function/(x)? 

Question (a) remains unanswered up to the present day. Question (b) is an¬ 
swered in the affirmative. The series (4) converges to /(x) if it converges at all. 
If we replace convergence in the usual sense by ( C 1) summability,then follow¬ 
ing Fejer, we can obtain elegant and satisfying results. Fortunately, for the ap¬ 
plications which we require, we can get by with very simple considerations. 

If 
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the infinite series (4) converges uniformly (and of course absolutely) for all x 
in (0,1) to/(x). 

It is clear that the convergence of the series in (5) implies the uniform and ab¬ 
solute convergence of the infinite series (4). Hence, we have 


S' 
J 0 


g(x y2™*dx = 


( 6 ) 


'k=-°o 

OO 

= w: 


,2uOac \-2 tx inx 


■) 


dx 


e 2nikx e ~2ninx ^ 


k=- oo 


(by virtue of uniform convergence) 


From this, it follows that f(x) and ^(x) are two continuous functions of x 
having the same Fourier coefficients; that is, 

00 f [f(x)-g(x)\e 27rinx dx = 0, 

J o 

for n =0, ± 1, • *. 

It is now plausible that the function f(x) - #(x) rnust be identically zero, and 
it is indeed true. Since this result can be established in a variety of ways, without 
further recourse to the theory of Fourier series, we shall accept it without 
furnishing any proof here. As usual, ingenuity is required if standard theory is 
not used. On the other hand, one of the basic results of the theory of Fourier 
series—the summability theorem of Fejer—yields the result immediately as a 
corollary of a far more general result. 

Exercise 

1. Consider the function + dy where s > 0 and 0 < a < 1. Ob¬ 

tain the Fourier series for a. (Hurwitz) 


3. The Poisson Summation Formula 

An important application of the theory of Fourier series is to the Poisson 
summation formula. This elegant and powerful technique can be used to derive 
a number of significant results in analysis. 




80 


Analytic Number Theory: An Introduction 6.3 


Let f(x) be a continuous function of x, defined for -°° < x < °°. Form the 
periodic function 

(1) £(*)=£/(* + ")• 

n = - oo 

For the moment, let us proceed quite formally, assuming that the series con¬ 
verges, and that the manipulations that follow are valid, and so on. Subsequently, 
we shall present a rigorous treatment. 

It is clear that 

(2) s00=^(x + l). 

Let us now invoke a mathematical principle first explicitly enunciated, and 
systematically exploited, by Hecke: A periodic function should always be ex¬ 
panded in a Fourier series. 

To obtain the Fourier coefficients of g(x) we write 



(3) = ^ J f(x + n)e~ 2lTlkx dx 

oo 

= y* f f(x)e' 2inkx dx =J f(x)e~ 2inkx dx. 

n=-°° 

Consequently, provided we can justify all of the above, we have the identity 

oo oo 

( 4 ) ^/(*+«)= ^ e 2mkx J f(x l )e~ 2ntkxi dx l . 

n=~ aa k~-°° 

The case x = 0 yields the Poisson summation formula: 

oo oo 

(5) £/(«) = 2 f^f^e^dx,. 

n = -°° k=-<*> 
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Exercise 

1. Evaluate 


f” («- [*])/'(*)<&> 

Jo 

use the Fourier series for (x - [jc] ), and thus obtain another proof of the Poisson 
summation formula. 

4. Some Simple Sufficient Conditions 

Examining the procedure of Sec. 3, we see that the imposition of some very 
simple conditions, which are easy to apply, will justify our procedures. Let us 
suppose that 

(a) the function f(x) is continuous for all real finite x; 

(b) the infinite series, 


£/(*+«) 

n—~°° 

converges uniformly in every finite x interval; 

(c) the infinite integral fZ\f(x)\dx converges; 

(d) the series 


2 i«*i. 

fc=-oo 

converges, where a k is the Fourier coefficient determined in 3.3. 

Then, under the foregoing conditions, the two sides of 3.4 exist and are equal 
for all x. This is Section 3 of this chapter. 


5. Application to the Theta Function 

Probably, the most famous application of the Poisson summation formula is 
to obtain the modular transformation. If we consider the sum £„ e"” 27It ,we see 
that it fulfills all the above conditions. Furthermore, the integral 
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can be evaluated, as we know. Consequently, we obtain the celebrated modular 
transformation. 

Exercises 

1. Apply the Poisson summation formula to the sum e ' n27Tt+2n7Tlz and thus 
obtain the transformation for the 0 function. 

2. Using the Laplace transform and the Hecke integral, derive from this for¬ 
mula the Iipshitz formula. 

3. From the Iipshitz formula, derive the functional equation for the 0 function. 

6. Multi-dimensional Version 

There is no difficulty in obtaining a multi-dimensional version of the Poisson 
summation formula. All that we need is the multi-dimensional analog of Fourier 
series. As before, we can give simple sufficient conditions for the validity of this 
formula. 

7. Interesting Domains 

In the application of the multi-dimensional formula it is important to have 
various domains of interest in analysis. In analytic number theory, we find the 
domain of Hecke, suggested by algebraic number theory, and the domain of 
Siegel, suggested by matrix theory. In analysis, we have the domain of Bochner, 

As we have shown in Chap. 1, other interesting domains arise naturally. 

8. The Evaluation of /ru e^ x>Ax ^dx 

For the application of the multi-dimensional Poisson formula to the multi¬ 
dimensional modular transformation, we must evaluate the /-<*> e~^ x,Ax ^dx. 
To accomplish this, it is convenient to use a little matrix analysis. As we indi¬ 
cate in the exercises, this is not necessary. 

We have 

A ~ T / Xl x 2 0 

0) I 

\ 0 X A 

where Xj • ■ • \ N are the characteristic roots. 
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Let us in the integral make the transformation x = Ty. Using the foregoing, 
the integral becomes 


( 2 ) *?+*' +K NyN 2 )dy, 

Using the one-dimensional case, we have the result 

(3) j ~ e - (y ’ Ay) dy = ^ 2 /\A |* 

where | A | represents the determinant of A. 

Exercises 

1. Show that the above formula is valid if A has the form B + iC where B is 
positive definite and C is symmetric and commutes with B. 

2. Show that the above formula is valid if B is positive definite and C is 
symmetric. 

3. Prove the formula in the text by reducing the quadratic form to a sum of 
squares. (Lagrange) 


9. The General Modular Transformation 

Using the multi-dimensional Poisson summation formula and the integral 
above, we can readily obtain the general modular transformation. 

What is difficult is in any particular case to pick out the theta functions that 
are pertinent. Here we are guided by algebraic number theory or by matrix 
theory in many cases. 

10. Sums of Squares 
We easily see that 



where r k (n ) is the number of representation of n as a sum of k squares. If we 
make the change of variable e~* = z, we can use Cauchy’s formula to express 
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r k (ri). We expect that the modular transformation will help us obtain an expres¬ 
sion for r k (n). To accomplish this, the Farey disection is very useful. 

This program was carried out by Hardy. In the general case, the same method 
can be used for Waring’s problem, as was done by Hardy and Littlewood. A dif¬ 
ferent procedure, based on trigonometric sums, was used by Vinogradov. 

11. Partition Functions 

Another important problem in number theory is the representation of an 
integer by integers of a given set. The most important example of this is the 
representation of an integer as a sum of integers. This is the famous partition 
problem, which has been studied extensively. 

A modular function again plays a central role, the modular function of 
Dedekind. 

12. Gaussian Sums 

One of the beautiful applications of the transformation formula of the theta 
function is to the evaluation of the Gauss sum, 


a) 


S(p,q)='£e-™ 

r=0 


where p and q are relatively prime integers. These sums are of great importance 
in number theory. 

To find the connection between this sum and the theta functions, let us take 
the function 


( 2 ) 



and examine its behavior in the immediate neighborhood of the line of con¬ 
vergence, Re(f) = 0. Set t = e + nip/q, where e is a small positive quantity and p 
and q are relatively prime positive integers. Then 


/(e + nipjq) = 1+2 

n =1 


(3) 
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= 1 +2 



2 

upon taking account of the periodicity of e~ mn p ^as a function of q. 
The function of e 


^ e ~(r+sq?e 
s=0 


e “W 2 e dw ^ 1 [~ e ~w 2 dw _ V* 

q q\Je Jo 2q\fe 

as e-*0. 

Hence, asymptotically, as e -+ 0, we have the equality 

( 6 ) 

2 

The periodicity of e 1tm p f q as a function of q enables us to write 

<7-1 <7 

(7) S(p,q)= 2 = 2' 

r=0 r~l 

We now employ the transformation formula 

(8) «')-(t)Xt) 

and repeat the process with the function f(n 2 /t). 

We have, for small e, 


(4) 

behaves like the integral 

(5) feWe, f 

J o J r 


n!_= 7,2 - A* ~ ™eIs} - e *? 2 + 0 ( c 2) 

t e + nip/q e 2 + w 2 p 2 /? 2 P 2 P 


Hence, as above, as e 0, 
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( 10 ) 



~ ~S(rq,p). 


Observe that the effect of the transformation formula of the theta function has 
been to invert the roles of p and q . 

Taking account of the fact that 



as e 0, we see that as e -* 0, 



Hence, equating the two asymptotic results. 


(13) 


_1_ 


qr-1 


. 2 , e ' ni/4 
~Ttir pjq _ e 


r =0 


\fp 


P -1 

1 

r=0 


nir 2 q/p 


a remarkable functional equation. 

Exercise 

1. Derive from the formula above the law of quadratic reciprocity. 


13. Eisenstein Series 


Consider the function defined in the upper plane. 


0 ) 



1 

(mz + n) k 


where the prime indicates that m and n are not both zero simultaneously. A 
series of this type is called an Eisenstein series. It is clear that/(z) satisfies the 
functional equation 

( 2 > m=z k f(iiz). 

More generally, we see that the function possesses a simple transformation 
rule under the transformation (az + b)/(cz + d). 
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Functions of this form arise naturally from the Poisson summation formula 
or the lipshitz identity. They also can be used to define general automorphic 
functions as was done by Poincar^. 

14. Functional Equations 

In the sections above, we put in the parameter additively. Now, let us put in 
the parameter in a multiplicative way. Consider the function defined by 


oo 

(1) F(x) =£/(«*). 

« = 1 

If we take the Mellin transform of both sides, we have 

( 2 ) M(F) = UsMf)- 

If we now invert, shift the contour, and use the functional equation for f (s) 
we obtain the Poisson summation formula. Naturally, the details require a bit of 
care. What is interesting about this method is that it shows that associated with 
every functional equation for a Dirichlet series is a summation formula. 

15. Use of the Laplace Transform 

We want a formula for the sum 

n 

(1) F(n) = £/(*). 

k=l 

We see that this function satisfies the equation 

(2) F(/i + l)-F(«)=/(« + l). 

Hence, we are led to consider the difference equation 

(3) F(x + 1)-F(x)=f(x), x>0. 

This difference equation can clearly be treated by use of the Laplace transform. 
We have 

e s ~l 


( 4 ) 
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If we employ the Iipshitz identity, and the inversion formula, we find again 
the Poisson summation formula. 

If we use the expansion of 


(5) 


1 1 
e s -1 s 


as a power series in s, we find a formula due to Fuler-Maclaurin-Abel-Plana. 

Exercises 


1. Use a similar method for the vector-matrix equation f(x + 1) = Af(x) + 
g(x). 

2. Use the same method for the equation u(x + 1) = u(x ) +/(*). 

3. Use the Mellin transform for the ^-difference equation u(qx ) = u(x ) +/(;c). 

4. Use the Lipshitz identity to find the nth derivative of l/(e* - 1). 

5. From this, obtain an estimate for the remainder term in the Taylor expan¬ 
sion of l/(e x - 1). 

6. From this, obtain an estimate for the remainder in the Euler-Maclaurin 
sum formula. 
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= i g t~ k/2 ' g ^P(n 1 , n2 , ■ ■ 

8 


,«/c)ex p 
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e 2irir 2 Jq _ 0 ~ l<? ) 

( 1-0 
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treatment of the sum Here d(n\ the Dirichlet divisor function, is 

the coefficient of n s in f 2 (S). 


n=l 


See 

Voronoi, C., Sur une Fonction Transcendante et ses Applications a la Summa¬ 
tion de Quelques Series, Ann. Sci. Ecole Norm. Sup. (Paris) 21 (3) (1904), 
207-267,459. 

For some interesting series connected with the Gauss circle function r(n) 
and the transformation formula 


t* ^ iiri)e~ nixt = f Vl 

n= 0 


5 ) r(n)e' nn/t , 

n=0 


see 

Kochliakov, S., Messeng. Math. 59 (1929), 1 . 

Further references will be found in the papers of Ferrar. A summary of recent 
results is given in 

Sklar, A., On Some Exact Formulae in Analytic Number Theory , Report of the 
Institute of the Theory of Numbers, Boulder, Colorado, University of Colo¬ 
rado, June 21 - July 17,1959. 

Detailed analysis is necessary for the sums that arise from these summation 
formulas. See 

Atkinson, F. V., The Mean Value of the Zeta-Function on the Critical Line , 
Proc. London Math. Soc. 47 (2) (1942), 174-200. 

Dixon, A. L. and W. L. Ferrar, Lattice-Point Summation Formulae , Quart. J. 
Math. (Oxford Series) 2 (1931), 31-54. 

Ferrar, W. L., Summation Formulae and Their Relation to Dirichlet y s Series , 
Compositio Math. 1 (1935), 344-360. 

-, Summation Formulae and Their Relation to Dirichlet *s Series. II, Composi¬ 
tio Math. 4 (1937), 394-405. 
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Bellman, R., An Analog of an Identity Due to Wilton, Duke Math. J. 16 (1949) 
539-545. 

See also, 

Pearson, T. L., Note on the Hardy-Landau Summation Formula, Canad. Math. 
Bull. 8 (1965), 717-720. 

Soni, K., Some Relations Associated with an Extension of Koshliakov’s For¬ 
mula, Proc. Amer. Math. Soc. 17 (1966), 543-551. 

Nasim, C., A Summation Formula Involving o k (n), k > 1 , Canad. J. Math. 
21 (1969), 951-964. 

Berndt, B. C., Arithmetical Identities and Hecke’s Functional Equation, Proc. 

Edinburgh Math. Soc. 16 (2) (1968/69), 221-226. 

Nasim, C., On the Summation Formula of Voronoi, Trans. Amer. Math. Soc. 
163 (1971), 35-45. 

Section 15 
See 

Chowla, S., On Some Formulae Resembling the Euler-Maclaurin Sum Formula, 
Norske Vid. Selsk. Forh. (Trondheim) 34 (1961), 107-109. 

For the Laplace transform see the book 

Bellman, R. and K. L. Cooke, Differential-Difference Equations, Academic 
Press, New York, 1963. 




7. Functional Equations 


1. Introduction 

In this chapter, we want to discuss the functional equations which determine 
the theta function. We have already seen how functional equations arise in the 
study of the gamma function and the zeta function. In the bibliography and 
comments at the end of the chapter, we will briefly discuss how functional 
equations arise in other parts of number theory. We shall give reference to work 
on differential difference equations. 

2. Functional Equations Satisfied by Theta Functions 

We begin with the function 


0 ) 


(z, o= 2 


e ~n nt^lnniz 


This function is obviously periodic of period one. This means that it satisfies 
the equation 

(2) 0(z + l,O = 0(z f t). 

To obtain a further equation, we use the obvious fact 

® 2*2 

n n +1 
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If we simplify the right-hand side, we obtain the functional equation 
(4) e(z + t,t) = e 2 ™e- 7tt O(z,t). 


Exercises 


1. Find the functional equation for O' /0. 

2. By using Rouche’s theorem, find the number of zeros of 0 in a fundamen¬ 
tal parallelogram. 

3. By evaluating f(Q'/0)e 2nniz dz around the fundamental parallelogram find 
the Fourier series for O'/d . (Hermite) 

4. Find analogous functional equations for 


2 2 

e ~(atn +2bmn+cn ) e 2mniz+2imiw 


~°°<m ,n<°° 


where the quadratic form is positive definite. 

5. Let /(/, g) be the Jacobian of the two functions / and#. Find the functional 
equations satisfied by /(/, g)/fg where / and g satisfy functional equations of 
the foregoing type. 

6. Find the Fourier expansion for /(/, g)/fg by considering the integral 


g)lfg)e^ 2minz+2nin ^dz dw, 


where the integration is over the fundamental parallelepiped. 

7. Show that the sum 


2 


_ 1 _ 

(m + in) 2k 9 


where the sum is over all values of m and n and the prime indicates that m and 
n are not simultaneously zero, can be evaluated using the theta function. 

8. More generally, evaluate the sum 


t 

S (m + nt) 2k ‘ 

(These sums are more easily evaluated using the Weierstrassian elliptic function.) 
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3. Determination of 0(z, t) 

What we now wish to show, conversely, is that these functional equations de¬ 
termine 0 (z, t). 

If we use the periodicity, we know that the function must have the form 


( 1 ) X U nV e2niriZ - 

n--°° 

If we use the second functional relation, we obtain a recurrence relation for 
u n (t). In this way we see that any solution of the two functional relations must 
be a multiple of 0 (z, t). The multiplying factor, however, may depend upon t. 

4. The Modular Transformation 

Consider the function 0 (z/t, l ft). To obtain the required transformation, we 
write 

(1) e 0 (, z/t,, 1 /t) = f(z, t). 

We see that the effect is to interchange the transformations z + 1 and z + t. 
Using the original functional equations, we see that /(z, t) satisfies the equations 

f(z + l,t)=f(z,t), 

( 2 ) 

f(z + t,t) = e' nt e 2mz f(z, t ). 

It follows that/is a multiple of 6. 

Exercise 

1. Write the zero of the theta function z = h(t). Show that h(t) satisfies the 
equation h(t ) = th(l/f). 

5. The Determination of the Multiplier 

It remains to determine the multiplier. As we shall see there are several ways of 
doing this. 

One way is to choose a particular value of z. For example, we can set z = 0. 
Thus, we reduce the general case to the special case. 
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6. Partial Differential Equations 

Another approach which is much more powerful in many cases, is to use the 
theory of partial differential equations. In particular, we use the theory of para¬ 
bolic differential equations, the heat equation or the diffusion equation. The 
fundamental observation is that 


2 

e ~n Titglnmz 

is a solution of the partial differential equation. 


(1) _ b 2 u 

~dt 3 

Since the equation is linear, a sum of two solutions is also a solution. If we use 
this fact on the left and on the right, we obtain the modular transformation. One 
advantage of this approach is that it can be used in many other cases. 


7. Generalizations 

As we have indicated in the exercises, the method can be extended to many 
other cases. Thus, we have a systematic way of handling theta functions of 
several variables. 

Also, the method can be extended to groups. We thus have a way of defining 
theta functions over groups. We shall not enter into this because it will take us 
too far from analytic number theory. 


8. fc-Dimensional Linear Spaces 

Consider the function /(z, t) which satisfies the functional equations 

{ l) f(z + t,t) = a(t)e 2k ™f(z, t). 

The same methods that were used above show that any function that satis¬ 
fies these relations is a linear combination of k particular solutions. 

If we consider the original functional equation, we see that the fcth power 
satisfies relations of the above form. We see then the reason why there must 
exist relations between particular solutions. 
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Exercises 

1. Obtain the functional equations for the square of 0 ( z , t). 

2. From this show that there must be a linear relation between any three 
solutions. 

3. How do we obtain the coefficients? 


9. The Case k-2 

Consider the four functions determined by taking all the combinations of plus 
or minus in the functional equations 

u(z + l,t) = ±u(z,t), 

u(z + t,t) = ± e ntt e 2mi u(z, t ). 

Let us prosaically call these functions f lt f 2 , fy, and / 4 . It will be recognized 
that fi is the function 6(z, t). 

It will be seen that the squares of these functions satisfy the equations 

, . f(z + \,t)=f{z t t\ 

( 2 ) „ . „ . 

f(z +1, t) = e‘ 27r, *e 4wiz n(z, t). 

Hence, there is a linear relation between any three squares. We can go further 
than this. Consider the functions f t (z + w, t)fj(z - w, r), where i and j runs 
through the numbers 1, 2, 3, and 4 independently. We see that these functions 
satisfy the equations in (2). Consequently, we have that this function is a linear 
combination of any two squares. The coefficients depend upon w. However, we 
see that this function is symmetric in z and w. Hence, we have 

(3) fi{z + w, t)fj(z - W, r) = afl(z, t)fl(w, 0 + bffo, t). 

The coefficients a and b depend upon i, j s k, and /. 

Exercises 

1. Consider the 16 functions determined by all combinations of plus or minus 
in the functional equations 

f(z + l, w, t x , t 2 , t 3 ) = ±f(z, w, t u t 2 , t 3 ), etc. 
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Prove that a linear relation exists between the squares of any five of these 16 
functions. 

2. More generally, consider the functions which satisfy 

g(z + l,t) = o/g(z, t), 

g(z + t,t) = cj s e mt e 2mz g(z, t), 

where co is a root of unity. 

3. By taking z = w, obtain a duplication formula. 

4. In general, obtain a formula for multiplication of the argument by n y 
where n is an integer. 

5. Obtain multi-dimensional analog. 

10. Complex Multiplication 

Consider the function d(zt, t). We see that this function has a very simple 
transformation under the transformation z + 1. Suppose we consider the trans¬ 
formation z + t and choose t to satisfy 

(1) t 2 =at + b, 

where the roots are complex, and a and b are integers. We see then that the func¬ 
tion d(zt, t) satisfies very simple transformations if t is chosen this way. This 
phenomenon is called complex multiplication. It was discovered by Abel and 
plays an important role in the theory of elliptic functions, algebraic number 
theory, and modular functions. 

Exercise 

1. Find analogs of complex multiplication for the higher-dimensional case. 

11. Doubly-Periodic Functions 

Consider the function defined by the quotient /}(z, r)//J?(z, r). Using the func¬ 
tional transformations, we see that this function is doubly-periodic. 

Thus, the theory of doubly-periodic functions may be made to depend upon 
the theta functions, as was done by Jacobi. One advantage of using the theta 
function is its rapid convergence. For theoretical purposes, in many ways the 
functions introduced by Weierstrass are superior. As usual, it is a good idea to 
have them both. 
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Miscellaneous Exercises 

1. What is the general solution of 

g(qd) = e—giq)! 

P 

2. Consider the infinite product f(q ) = n„(l + qx n ). Show that f(q) satisfies 
the functional equation f(qx) = (1 + x)f(q ), and thus determine the power 
series expansion. 

3. Consider in a like fashion the infinite product f(q) = n„(l - qx n ). 

4. Obtain a functional equation for the infinite continued fraction 


u(q) = 


X 

1 +qx 


1 + q 2 x 


5. Show that this nonlinear qr-difference equation can be converted into two 
linear qr-difference equations by means of the transformation u(q) = v(q)/w(q). 

6 . Let /(m, n ) be the number of steps in the Euclidean algorithm where 
m > n > 1. Prove that f(m, n) satisfies 

f(m, n) = 1 +/(«, m - [m/n] n). 

See, 

J. D. Dixon, A Simple Estimate for the Number of Steps in the Euclidean Al¬ 
gorithm , Amer. Math. Monthly 78 (1971), 374-376. 

7. Consider the Newton interpolation formula 

f(n) = a 0 + a v n + a 2 - + * * *. 

Show that the coefficients may be determined by the operation 

8 . Show that a continuous limit of this formula is the Taylor series. 

9. Consider the analog of the foregoing interpolation series 

f(p) = a o + a 1 7~TZ + '• 

(«-l) 
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Show that the coefficients may be determined by the operation 

;/(„). 

10. Use the foregoing interpolation formula on the function n"J^(l - q k ). 
(See, 

Bellman, R., A q-Version of the Newton Interpolation Formula and Some 
Eulerian Identities, Boll. Un. Mat. Ital. 16 1961,285-287.) 

11. Consider the function f(x) = 2 £ =0 x 2 ” . Prove that this function satisfies 
the equation f(x) =x +/(xr 2 ). 

12. Prove that this functional equation determines the power series. 

13. Consider the sequence 

= a n + b n 
" +1 2 ’ 

^n+1 — > 

a 0 = l, b 0 =x, a<x <<*>. 

Prove that we have a n > b n . Show that the sequence a n is monotone de¬ 
creasing, and the second sequence b n is monotone increasing. Hence, show that 
they both converge, and that they converge to a common limit 0 (x), the arith¬ 
metical geometric mean of Gauss. 

(This mean plays an important role in the evaluation of elliptic integrals, the 
Landen transformation. See 

Whittaker, E. T. and G. M. Watson, A Course of Modern Analysis, 4th ed., 
Cambridge University Press, Cambridge, 1935.) 

14. Show that <p(x) satisfies a modular transformation. 

(This function can be made the basis of a theory of elliptic functions. See 

Gauss, C. F„ Gesammelte Werke. 

15. Prove that 

a«+i -b n+ 1 =(vs;-v^) 2 . 

Hence, show that we have quadratic convergence. 
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Rogers, L. J., Proc. London Math. Soc. (II), 16 (1917), 315. 
and for another class of related functions, see 

Basoco, M. A., On Certain Arithmetical Functions Due to G. Humbert , J. Math. 
Pures Appl. 9 (1947), p. 237. 

The method of functional equations enters into many parts of analytic num¬ 
ber theory. In particular, differential difference equations occur. 

Differential difference equations play a prominent role in number theory, 
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showed by means of a Landen transformation that 

/ V_ do _ = r*/2_ dO _ 

o ( a\ cos 2 0 + b\ sin 2 0) % J o ( a 2 cos 2 0 + Z? 2 sin 2 0) % ’ 



7.9 Functional Equations 


105 


a result which yields, upon repeated application, the relation 

r*l 2 _ dd _ = tt/2 

J o (a 2 cos 2 0 + b 2 sin 2 0) % M(a,b) 
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2 2 2 
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manifestations of a general principle that partial differential equations have 
similarity solutions for small times and separation-of-variables solutions for 
large t. 

The functional equation for the theta functions, particular solutions of the 
heat equation, is a special case of a far more general transformation formula 
valid for general solutions of the heat equation. 

Let u(x , y, z, t) be a particular solution of the three-dimensional heat equation 


1 du _ d 2 u + b 2 u + d 2 u 
?d7 dx 2 by 2 te 2 ' 


Then the function 


U(x, y, z, t ) = f 3/2 u 



is also a solution. There are corresponding results for the ^-dimensional heat 
equation, and analogous results for the wave equation. 

We have presented the foregoing method of establishing transformation for¬ 
mulae, since it has been used by Maass and Siegel to establish some important 
identities which to date have not been derived by any other means. See 

H. Maass, Uber eine neue Art von Nichtanalytischen Automorphen Funktionen 
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Klasse, Berlin, 1927, pp. 158-161. 
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Section 10 

Complex multiplication played an important part in the “Jugendtraum” of 
Kronecker. 

Section 11 

An excellent discussion of the doubly periodic function is in the book by 
E. T. Whittaker and G. M. Watson. See 

E. T. Whittaker and G. M. Watson, A Course of Modem Analysis , 4th ed., 
Cambridge University Press, Cambridge, 1935. 

For the theory of multiply periodic functions, see the book 

G. F. Baker, An Introduction to the Theory of Multiply Periodic Functions , 
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8. The Euler 0 Function 


1. Introduction 

In this chapter, we will study one of the fundamental functions of number 
theory, the Euler \p function. 

First, we will establish the fundamental multiplicative property. Using this, we 
obtain an Euler product from which we will derive a very simple representation. 
Using this representation, it is easy to derive various mean values. 

Finally, at the end of the chapter, we shall discuss the Ramanujan function. 

2. The Multiplicative Property 

In this section, we shall derive the fundamental multiplicative property of the 
Euler <p function. We use a simple but useful technique which appears in many 
parts of number theory. We calculate our quantity by two different methods and 
equate the results. 

It is convenient to calculate the number of integers, <mn, which have a fac¬ 
tor in common with the product mn (where m and n are relatively prime). 

On one hand, this number is 

(1) mn - ip(mn). 

On the other hand, this number is 

(2) (m - v(m))n + m(n - p(n)) - (m - <p(m)) (n - y>(n)). 

The third term arises because the first two terms count several quantities 
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twice. Equating the two expressions we obtain, after some simplification, the 
fundamental result 

(3) v(mn) = ip(m)ifi(n), (m, n) = 1. 

It is clear that 

(4) 'P(p a ) = P a ~ 1 (p- 1). 

Hence, using the fundamental theorem of arithmetic, we have 


(5) 

p\n 

Hence, we have 


( 6 ) 

Exercises 



1. Show that lim„-+oo <p(n)/w = 1. 

2. Show that lim^oo ip(n)/n = 0. 

3. Show that if (m, n) = 1, 1 < m < n, then n - m is a number of the same 
form. 

4. Using this fact, evaluate the sum £( m , n )=i m where 1 < m < n. 

5. Define S k = S( m tH y=i m k , where 1 < m < n. By considering the sums S 3 
and *S 4 evaluate S 2 and S 3 . 

6 . Generalize, to evaluate S k . 

7. Show that 


n=l n k =1 


I4Q x k 

k ( l-X k )’ 


3. The Euler Product 

Using the fundamental multiplicative property, we have the relation 
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( 1 ) 



+ 


+ 



) 


■n( i+ s= a / i - i " , ' ) ) 



= r(s_j ) fori?e(s)>2 

f(s) 


From this we derive the representation 





The quantity <^(«)/« is so useful to consider that we shall call it f(n). 
Exercise 


1. Show that the two representations of ip(ri) are equivalent. 
4. The Mean Value of ip(n) 

Using the representation above, we have 

o) 5/(«) = 5YX • 

n -1 w=l\/c|« / 

Inverting the order of summation we have 
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using the relation 

(3) 0 - [x] < 1, 

we have 



This is a remarkably sharp estimate. Using the estimates we derived in Chap. 1, 
we can write this as 

(5) ^ + 0(logA0. 

7r 

Exercises 

1. Derive an estimate for y(k). 

2. Derive the constant represented by the symbol 0. 

5. The Mean Value of <p(/z) 2 

To obtain the order of magnitude of the mean value of f(n ) 2 , we proceed as 
follows. We write, 


( 1 ) 


N N / 

^m 2 = ^m(YlJi(r)/r\. 

k =1 k=l \ r\k / 


Inverting the summations, this may be written 

N 


( 2 ) 




k=l 


n^O(k) 
<1 <k<N 


Using the results of the previous section we see that the required mean value 
has the form 


(3) 


c 1 iV + 0(logA' 2 ) 


where is a constant. 
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Exercises 

1. Evaluate in this way the sum if(k ) 3 . 

2. Find the mean value of ip(n) 2 . 

3. Find the order of magnitude of the Ingham sum /(«)/(« + a). 

4. Find the order of magnitude of the Ingham sum 2% =1 <p(n)ip(n + a). 

6. An Alternate Method 

It is clear that we may proceed inductively to find the order of magnitude of 
the sums 

N 

(i) 2 mk - 

n-\ 


However, the arithmetic soon gets rather unpleasant. Let us then give a dif¬ 
ferent approach. 

We have 



= ?(s)s(s), 


where g(s) has an absolutely convergent Dirichlet series for s > 
Hence, we may write 


/(«)* = £ v 

r\n 


( 3 ) 
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where a r is absolutely convergent. Thus we have 

( 4) £/w*=i4fl~*i> 

n =1 r= 1 L J r=l 

We leave it as an exercise to evaluate the error term. 

Exercise 

1. Find the order of magnitude of £^ =1 ip(ri ) k . 

7. The Mean Value of 

In what follows, we assume that p(n) has no multiple factors. 

Let us now consider the mean value of ip(p(n)) where p(n) is a polynomial. 
We have 

N N / \ 

(1) y /ooo)=y (y mw/*). 

n=l n=l fclp(n) 

Inverting the order of the summations, we have 

if s 1 

fc=l /(«)=0(fc) 

l<n<iV 

(2) -£«S.rfl 40008^ 

k=l L J 


8. The Ramanujan Function 

Let us now consider an important function which may be considered a gen¬ 
eralization of p(n). 
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We define 

(1) c q (n)= X 

(m,qr)=l 


2nimn 
e - 

Q 


This function has some interesting properties which we shall discuss further 
in Chap. 9. 
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Miscellaneous Exercises 

1. Prove that c q (n) is a multiplicative function. 

2. Evaluate c q (n) in terms of <p(n). (Holder) 

3. Prove that the c q (n ) are orthogonal. (Carmichael) 
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9. The Divisor Function 


1. Introduction 

In this chapter, we study the divisor function. 

In the second section, we recall some facts about this function. In the third 
section, we consider its mean value. In the fourth section, using a simple geo¬ 
metric method, we show that the error term can be considerably improved. In 
the fifth section, we show an analytic equivalent of the geometric result which 
we use in a section below. In the sixth section, we consider the mean value of 
d(n) 2 . In the next section, we show the connection with certain Diophantine 
equations. In the next section, we give the famous Perron sum-formula. Then we 
show how this formula may be extended by a summability technique. Then we 
turn to the mean value of d(an 2 + bn + c ). In the next section, we show why 
the method does not work for d(p(n )), where p(n) is a general polynomial. In 
the final section we show how Ramanujan expansions can be used for the 
evaluation of the order of magnitude of o(p(n)). 

2. The Divisor Function 

Let us now consider some facts about the divisor function. This function is 
the number of the divisors of n . Analytically, we may write 

(1) d(n) = ^l. 

k\n 


A generating function is 
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(2) rto-i*?. 

n =1 

Since d(n) is multiplicative, we have 


(3) 

4«)=[]( 1 +a <). 


p 

where 


(4) 

n = npf«. 


In this chapter we shall not use the multiplicative result. 

Exercises 


1. Show that d(ri) = 0(« e ) for any e > 0. 

2. Show that there are integers n such that d(n) > (log ri) k for any k greater 
than 0. 

3. Find the order of magnitude of d(ari). 

4. Find the order of magnitude of 2^ =1 d(an + b ). 

5. Find the order of magnitude of 2^ =1 r(n). 

6. Show that 


^ x n 

^ n s ( 1 - x n ) 
n= 1 


oo 

*2 a s( w >*"- 


n=l 


7. Show that o s (ri) = n s o. s (n). 


3. The Mean Value of d(n) 


We have 


(i) 


N N 

5 >- 2 ( 2 ‘) 


« = 1 «=1 fcjrt 

Inverting the orders of summation, we have 


iv 




M — 1 



( 2 ) 
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a deceptively simple relation. 
Thus we have 


(3) 


N' N 

£d(«)=;v£l+0(JV). 

fl=l *=1 


Using the fundamental relation of Euler, we have 


(4) 

I 


N 

^ d(n)=NlogN + 0(N). 

n=l 


Exercises 

1. Find the order of magnitude of 2^ =1 d(an). 

2. Find the order of magnitude of d(an + b ). 

3. Find the order of magnitude of o s (n). 

4. Show that k and n/k are simultaneously divisors of n. These are distinct if 
n is not a perfect square. 

5. Using this fact, obtain the order of magnitude of a 5 («) with an error 
term 0(y/n). 

4. A Simple Geometric Approach 

We can do much better. 

Consider the following diagram 



i 

i 


I 
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We know that the number of the lattice points determined by 1 < xy < n 
is given by the expression [n/k]. From considerations of symmetry 

the number of lattice points in the region I is the same as in the region III. 
Hence, we have the relation 


( 1 ) 


« r- -i V” 
k=l k=l 


Hence, we have 

(2) ^ d(n) = N\ogN+ (2c - 1 )N + 0(\//V). 


N 


n=1 


Exercises 


1. Find the order of magnitude of (x/n - [x//i]). 

2. Use the same method to find an error term for E^ =1 r(n ). 

5. Analytic Equivalent 

An analytic equivalent of the simpler geometric method presented above is the 
observation that k and n/k are simultaneously divisors of n. Hence, we have the 
relation, which we will use below, 


( 1 ) 


where 


( 2 ) 


d(ri) = 2 1 - e(n), 

k\n 

k<sfn 


e(n) = 1, if n is a perfect square, 
= 0, otherwise. 


6. The Dirichlet Divisor Problem 

Let us write 

N 


A (AO = ^ d(ri) -yVlogjV- (2c - 1)N, 


n =1 


(i) 
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Here Cis the Euler Constant. 

The Dirichlet divisor problem is the estimation of A(w). This is very difficult. 
Two methods are used at present. The first is due to Voronoi (Woronow). The 
second method uses trigonometric sums. Both involve complicated analysis. 

In Chap. 12, we show that Tauberian theorems can be used to yield order of 
magnitude for the mean value of A (n ) 2 . 

The same problem for the sum 2^ =1 d k (n) is called the Piltz divisor problem. 
Here d k (n) is a higher divisor function. 

7. The Mean Value of din) 1 


We would expect from the foregoing that d(n ) behaves like log n. To our sur¬ 
prise, this is not so. This means that occasionally, as we have pointed out d(n) 
has very large values. To show this, let us find the mean value of d(ri ) 2 . 

We have 


( 1 ) 


2 d(w)2= 2 ^(S 1 )- 

n=l n -1 k\n 


Inverting the order of summation, and evaluating the sums that arise, we find 
that 

N 

(2) ^ d(n) 2 = CjTVCiogAO 3 + 0(N(logN ) 2 ). 

n -1 


We leave the details to the reader. 

Exercises 

1. Use the generating function to obtain the foregoing results. 

2. Evalute 2^ =1 d(n ) 3 . 

3. Generally, evaluate 2^ =1 d(n) a , where a is a positive integer. What happens 
if a is not? 

4. Find the order of magnitude of 1^=1 o s (n)o s (n + a). 

5. Evaluate the constant C x by two methods and equate. 

8. Connection with Diophantine Equations 

We observe that d(n) is the number of solutions of x x x 2 = n. Consequently, 
the number of solutions of 
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(1) *l*2 + *3*4=« 

may be written as 


N -1 

(2) ^jT d(k)d(n - k). 

k=l 


Similarly, the number of solutions of x x x 2 - x 3 x 4 = n, for x x x 2 is equal 
to; 

N-l 

(3) ^ d(k)d(n + fc). 

Ar=l 


We call sums of this type Ingham sums. 

Exercises 

1. Find the order of magnitude of the Ingham sum S^ =1 d(n)d(n + a). 

2. Find the order of magnitude of d(n)d(N - n). 

3. Find the order of magnitude of r(n)r(N - n). 

4. Find the order of magnitude of r{ri)r(n + a). 

5. Consider the number of solutions of x\ + x\ +xl + x\ =n. Show that this 
number may be written as S^= 0 r(n)r(N - ri). 

6. Evaluate the order of this sum, and hence prove Lagrange’s theorem that 
every integer is the sum of four squares. 

9. The Perron Sum Formula 


Let us derive an expression for the partial sum of the coefficients of the 
Dirichlet series. This expression will be very useful, as we shall see in deriving 
various mean values. We start with the integral 


( 1 ) 



/ 


x s ds 

s 


Here Cis a line parallel to the imaginary axis and lying to the right of it. We have 



9.8 The Divisor Function 


/(*)= 0, 0 <X < 1, 

= 1, X>1. 
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( 2 ) 


The first case is obtained by shifting C to the right; the second case is obtained 
by shifting C to the left and taking account of the residue at 5 = 0. Notice that 
we say nothing about the value for x = 1. 

Using this integral, we have formally the result, where we assume x is not an 
integer, 

(3) Y = ~ f~mds, 

^ 2m c s 

l<n<x 

Here C is a straight line to the right of the singularities of f(s). Where 

oo 

(4) 

n=l 

It remains to prove this formula. We write 

N 

(5) /(*) = 2^ +7,Jv(s) 

n-\ 


where N is larger than jc. Thus, we have 


( 6 ) 


_L f fcW ± - l fx 1 y 

2m Jc s 2mJc s +* r? 

n -1 

. y 

^ 2m Jc 


1 <n<x 


+ J_ f x s r N (s)ds 
2iri Jc s 

x s rfif(s)ds 

s 


It remains to estimate the second term. To do this, we shall integrate by parts. 
We have 


( 7 ) 


fx s r N (s)ds= 2 
n=JV+l 


log (n/x) 
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We may integrate term-by-term if C is far enough to the right. If N is large 
enough, the term log (n/x) causes no difficulty. Thus, we see upon carrying out 
the integration by parts, that we can make the integral as small as desired by 
taking jV sufficiently large. Thus, we obtain (1), the Perron sum formula. 


Exercises 


1. From the representation 


x<n<x+y 




S 2 (s)ds, 


can we find a better estimate for the sum than the estimate given above? 

2. Shift the contour, use the functional equation for f (s), and thus derive 
the representation of Voronoi for the error term. 


3. Do the same for r(n). 

4. Do the same for r k (n). 

5. Use the Perron sum formula for the interval (x + y,x). By choosing y as a 
function of x, we can obtain a better result than by using the standard error 
term. 


6. Show that the Perron sum formula can be obtained by taking the Mellin 
transform of the function 2> n <x 

7. What do we get if we apply the Laplace transform? 

8. Obtain a similar result for m<x, n<y. 


10. Logarithmic Summability 

It is easy to derive a generalization of this formula. 
We begin with the representation 

( 1 ) (logx)*' 1 _ 1 r x s ds 

1 } (k- 1 )! 2m Jc s k 


for any positive integer k. This can be proved in exactly the same way as the 
formula above. Using this formula, we have 



9.11 The Divisor Function 


133 


5) 0°g x/n) k a„ 

tl<X 


kj_ 

2in 


r * s m 

Jc s k+l 


ds. 


Exercises 

1. Use this result to obtain an error term for the divisor function. 

2. Use this result to obtain an error term for the circle function. 

3. Use this result to obtain an error term for the higher divisor function. 


11. Asymptotic Behavior 

The Perron sum formula is used in the following way. We shift the contour to 
the left and determine the residue of the integrand. If we know the analytic be¬ 
havior of f(s), this is simple. This residue gives the order of magnitude of the 
partial sum. It remains to estimate the integral term. As we shall see, the mean 
value of the function f(s ) can be used for this. In general, we do not obtain a 
very sharp estimate in this way. However, we get the order of magnitude very 
easily and in Chap. 12; we shall see that this gives us an estimate for the mean 
value of the square of the error term. 

Let us now give an example of how the sum formula can be used. Consider 
the expression. 


a) 


y d(n>-±- 

^ 2m Jc s 

1 <n<x 


Here C is a line lying to the right of the line a = 1. Shifting the contour to the 
left, we pass a single pole at s = 1. 

Thus we have the expression 


( 2 ) 


A(x) = 


2m 



Here C' is a line to the left of a = 1. In order to make this integral as small 
as possible, we would like to shift as far as possible to the left. To estimate 
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this remainder term, it is convenient to use a contour of the following form 



It remains to estimate the integral along the various parts. The integral over the 
lines parallel to the jc-axis are easily estimated. The integral along the line to the 
left of the residue at s = 1 is estimated by using a mean value result for f (s). 

We have 

s I •'Is 

The integral here is taken over this vertical line. Using the Cauchy-Schwarz 
inequality we have that this term is bounded by the expression 

If I / V \s\) 




Thus, we see that our ability to estimate the error term depends upon our 
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ability to estimate 


(5) J l?(s)l 4 <fc. 

Thus we see how the estimate of the error term can be made to depend on the 
estimation for various mean values. In Chap. 12, we will pursue this subject 
further. 

12. The Mean Value of d(an 2 +bn+c) 

Using the representation of d(n) given above, in (1) of Sec. 5, we have 

1 

+bn+c) 

\an 2 +bn+c\ 

Inverting the order of summation and using some algebraic number theory, we 
have the estimate. 


a) 


N N 

5) + bn + c) = 2 


n=l 


n=l\k/(an z 

k< 


N 

(2) ^ d(an 2 + bn + c) = f(a,b f c ) N(logN ) + 0(AT) 

n =1 

Exercises 

1. Find the order of magnitude of 2^ d(n 2 ). 

2. Find the order of magnitude of o s (an 2 + bn + c). 

13. The Mean Value of d(p(n)) 

We assume, as usual, that p(n) is irreducible. 

If we try the same method for the mean value for d(p(n)) > we find that the 
error estimate overwhelms what is supposed to be the principle term. We conjec¬ 
ture that 

N 

(1) ^ d(p(n)) = c 2 N\ogN + 0(AO, 

«=i 

where the constant c 2 depends upon the polynomial pin). 


136 


Analytic Number Theory: An Introduction 9.10 


14. Ramanujan Expansions 

If we wish to obtain more precise results, it appears that we must use a differ¬ 
ent method. Fortunately, we have expansions which permit us to obtain more 
refined mean value estimates. Such expansions were given by Ramanujan. Thus, 
we have, 


0 ) 


0-s(«) = t(l + «)]£ 

< 7=1 




Re(s) > 0, where, as usual o~ s (ri) = X k \ n l/k s . 

Using some estimates for trigonometric sums various mean values can be ob¬ 
tained from this expression. Reference is given at the end of the chapter. 

Exercises 


1. Obtain the foregoing expansion by decomposing the function 


oo 

2 


n s { 1 -x n ) 


into partial fractions. 

2. Obtain an expansion for II^ =1 (1 +x”). 
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Miscellaneous Exercises 

1. Show that every number is uniquely representable by powers of two. 

2. Show that an analytic equivalent of this statement is the identity 

— = (1 +x)(l +* 2 )(1 +x 4 )---. 

1 -X 

3. Establish this identity by means of the identity 

2 n 

— —— = (1 + x) (1 +JC 2 ) • • • (1 +X 2iV 1 ). 

1 -X 

4. Let a(n) be the number of ones in the representation of n. Find the order 
of magnitude of ct(ri). 

5. Find the order of magnitude of 2^ =1 a(a(«)). (See 

Bellman, R., and H. N. Shapiro, On a Problem in Additive Number Theory, 
Ann. Math. 49 (1948), 333-340. 

Grosswald, E., Properties of Some Arithmetical Function, J. Math. Anal. Appl. 
28 (1969), 405-430. 

Stewart, B. M., Sums of Functions of Digits, Canad. J. Math. 12 (1960), 374- 
389. 

Mahler, K., On the Generating Functions of Integers with a Missing Digit, K’o 
Hsiieh (Science) 29 (1947), 265-267. 

Alexander, R., Remarks About the Digits of Integers, J. Austral. Math. Soc. 
12(1971), 239-241. 

Thorp, E., and R. Whitley, Poincare’s Conjecture and the Distribution of Digits 
in Tables, Compositio Math. 23 (1971), 233-250. 

6. Can one find the order of magnitude of d(d(nj)l 

7. Show that 


Y_!_ 

4* d(ny 
n-l 


= f(s)V), 


where g(s) has an absolutely convergent Dirichlet series for RE(s) > 
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10. The Squarefree Problem 


1. Introduction 

In this chapter, we want to study some aspects of squarefree numbers. 
Squarefree numbers are numbers which possess no repeated prime factors. 
We would like to solve the corresponding problems for primes. What we hope 
is that consideration of squarefree numbers will give us some ideas of the order 
of difficulties of the corresponding problems for primes. What we find is that 
even the problems for squarefree numbers are quite difficult. 

In the second section, we find a representation which makes the mean value of 
squarefree numbers quite easy to obtain. When we try to extend this method 
to numbers of the form n 2 + 1, we find, as usual, that the estimate of the error 
term causes difficulty. We can circumvent this difficulty by considering the Pell 
equation, an interesting equation in its own right. Then we turn to polynomials 
of higher degree. In the general case, all we can do is make a conjecture. The 
method used here is too weak. What we would like to prove is that there are in¬ 
finitely many squarefree numbers of the form n 2 + 1 without finding a mean 
value. 

We would like to consider the general quadratic polynomial an 2 + bn + c. 
We consider the polynomial n 2 + 1 in order that we can use the known prop¬ 
erties of the Pell equation. This simplifies the proof. 

2. The Generating Function 

We have « 

n~ 1 p 
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P 

= f(s)/f(2a). 

From this, we obtain the fundamental representation 

(2) M 2 (») = 2 V( k )- 

k 2 \n 

Exercises 

1. Obtain this representation by argument without using a generating function. 

2. Show that the generating function of the rth power free integers is given 
by f (s)/f (rs). 

3. The Mean Value of n 2 (n) 

We have 


( 1 ) 


Exercises 


*■<*)) 

n =1 n=1 fc^ln 

') 

k=l n=0(k) 

l<n<N 



K 

k=l 

= ^ + 0(VAT) 


1. Obtain a corresponding result for rth power free integers. 
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2. Improve the error estimate by assumption of the Riemann hypothesis. 

3. From the representation 


rt . 2m Jc\ 1 ' i(2s) 


x <n<x+y 


can we find a better estimate for the sum than using the estimate of the error 
term above? 


4. The Mean Value of ii 2 (n 2 + 1) 

If we use the same technique as above, we find we have difficulty, as usual, 
with the error estimate. Consequently, we must employ a new technique. 

We shall study the Pell equation to obtain an estimate for the number of 
solutions of 


(1) n 2 +1 = 0(> 2 ) 

for n < x and y sufficiently large. 


5. The Pell Equation 

The congruence above is equivalent to the equation 

(1) « 2 + l=fcy 2 . 

The equation 

(2) n 2 - km 2 = -1 

is known as the Pell equation, a very interesting equation in its own right. Let us 
study some of the properties of the solutions of this equation. 

We associate with every solution the number n + my/k. We have the basic 
factorization 

(3) n 2 - km 2 ~ (n -my/ic) (n + my/k), 

which shows the genesis of the association above. 

It follows from this factorization that the product of two solution numbers 
yields another solution number. Even further, the quotient of two solution 
numbers yields another solution number. What we want to show is that all 
solutions are generated by one solution. 
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To do this, we proceed as follows. Associate, as we have said before, with the 
solution the number n + m y/k. If n and m are positive this quantity is greater 
than one. Let n x + m x y/k be the smallest positive solution greater than one. 
Suppose that there is another solution n 2 + m 2 y/k which is not a power of the 
first. Then it must be between two powers of the first. 

By this we mean 

(4) (n +m y/k) a <n 2 +m 2 y/k < (n x + m x y/k) a + l . 

Consequently, we have 


(5) 


j < n 2 +m 2 y/k 
(n x + m x y/k) a 


<n t 


+ m 1 y/k. 


From what we have said before, the quantity (n 2 + m 2 y/k)/(n x + m x \/k) a 
represents a solution of the Pell equation. On the other hand we see from Eq. 
(5) that it lies between one and n x + m x y/k , a contradiction. Consequently, all 
solutions of the Pell equation are powers of the solution equivalent to n x + 
m x y/k. 


Exercises 


1. Show that all the numbers equivalent to the solution of the Pell equation 
form a group with a single generator. 

2. Associate with the solutions of x 3 + ky 3 + k 2 z 3 - 3kxyz = 1, the number 
x + $k y + $k*z. Show that these numbers form a group. 

3. Does this group have a single generator? 

6. The Number of Squarefree Numbers of the Form n 2 + 1 

We are now ready to discuss the number of squarefree numbers of the form 
n 2 + 1. The congruence we obtain leads to the equation n 2 + 1 = ky 2 . We could, 
if we wanted to, use the foregoing methods to discuss this equation. Instead, 
let us use the fact that the quotient of two solutions of this equation is a solu¬ 
tion of the Pell equation. If we consider the interval 1 <«<:*:, 1 
we see that for fixed k the number of solutions in this interval is 0(logx). Con¬ 
sequently, to estimate the sum 2^ =1 p 2 (w 2 + 1), we proceed as follows. We 
write 

^ V (« 2 + 0=2 ( 2 m (* o ). 

n= 1 n=l fc 2 |(n 2 +l) 


0) 
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Inverting the order of summation, we break up the resulting sum into two 
parts. Thus, we have 


N 2 N N 2 

® 2 - 2 * 5 , • 

n=l k =1 k^N^+l 

In the first sum we proceed as before. To estimate the second sum we use 
the fact that the number of solutions of n 2 + 1 = ky 2 in the interval 1 <n<N, 
1 < k < N is 0 (log N) for each fixed y. Consequently, if y is in the interval k 
> N 2 & + 1, we see that k must be in the interval 1 < y < N. We see that the 
sum is at most 0(N 2 P\ogN). Consequently, we have the result 

N 

(3) 2 M 2 (« 2 + 1) = CyN + 0(N 2,3 logN). 

n =1 

The equation need not have a solution for every k. Obviously, it is not neces¬ 
sary to consider these equations. 

A very simple sieving argument due to H. N. Shapiro works for the case of n 2 
+ 1. Namely: let x be an integer > 4, Q(x) = number of n 2 + 1, n = 1,2,..., 
x which are squarefree, B(p 2 ) = the number of solutions modp 2 to n 2 + 1 = 
0 (mod p 2 ), p a prime. Then we know that .5(2 2 ) = 0, and for p an odd prime, 
B(p 2 ) = 0 for p = -1 (mod4) and B(p 2 ) = 2 for p = 1 (mod4). Then 


(4) Q(x) >x - ^ (no. of n <x s.t. n 2 + 1=0 (modp 2 )) 

P< X 

so that 

(5) Q(x)>x- J ^B(p 2 )-2 J • 

p< \[x sfx <p<x 

But using the information on B(p 2 ) we have 



and from elementary prime number theory 2 1 < c(x/ log*), so that (2) 

yields 
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(7) g( x )>*A _£_'\>r!L (for jc large), 

\2 log*/ 3 

and this certainly implies infinitely many squarefrees of the form n 2 + 1. 

This simple argument can be elaborated upon so as to give an asymptotic 
result. Namely, let L = L(x ) to be chosen later, and let A = n p< x p. Then 


Q(x ) = (number of n <*x such that n 2 + 1 is not divisible by any 



k 2 /n 2 +1 



Then it is easy to see that 


so that putting everything together 



Thus taking L=c x (log*) (loglog*), (for small c x > 0), we get 
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Exercise 

1. Find the order of magnitude of Z^ =1 fi 2 (an 2 +bn+ c). 

7. The Number of Squarefree Numbers in a General 
Polynomial Sequence 

If we try to employ the foregoing method for the general case, we encounter 
difficulties. Consider, for example, the case where we have n 3 + 2. We know, 
from Thue-Siegel-Roth theorem, that the equation n 3 + 2 = ky 2 has only a finite 
number of solutions. However, at present, we have little information about what 
values of k have a solution at all, and how to estimate the number of solutions 
in this case. The method of Roth can be used to estimate the number of solu¬ 
tions but this estimate is too weak to use. 

Using difficult techniques, particular polynomials can be treated. 

It is possible that a method based on trigonometric sums will provide the in¬ 
formation we require. In any case, we conjecture that if p(n) is irreducible, we 
have 

N 

(1) ^M 2 (p(«)) = c 2 ^ + 0(AO, 

n =1 


where c 2 is a constant which depends upon pin). 

The correct condition is not that p(n) be irreducible but that there does not 
exist a prime q such that q 2 divides p(n) for all n. Then one conjectures (1) 
with a positive constant c 2 . (Note that an irreducible polynomial can be divisible 
by a square factor for all integer values of argument; e.g., x 4 -x 2 + 4 is always 
divisible by 4.) 

If f{x) = II/ =1 gfa) where each^(jc) is either a linear or quadratic polynomial, 
then for x - 1, • • • , N the number of squarefree fix) is cx + 0(xr) for some 
c > 0. 




150 


Analytic Number Theory: An Introduction 10.6 


Miscellaneous Exercises 

1. Find the density of integers n where n and n + 2 are both squarefree. 

2. Find the density of squarefree numbers in an arithmetic progression. (Landau) 
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11. The Prime Divisor Function, Selberg's 
Sieve Method, and Algebraic 
Independence 


1. Introduction 

In this chapter, we shall consider three different topics, the prime divisor 
function, Selberg s sieve method and algebraic independence. We shall first ob¬ 
tain the Euler product for 2^ =1 2 w ^/n s , where oo(n) is the number of prime 
divisors of From this, it is easy to obtain a representation which allows us 
to derive the mean value. Then, using the fact that c o(n) is an additive function, 
we shall obtain its mean value. Then we obtain the mean value of the square of 
co(ri). With the aid of these results we can get a meaning to the statement about 
the average number of prime divisors of an integer. 

Next, we turn to the sieve method of Selberg, giving his original proof. This is 
one of the most elegant and powerful methods in analytic number theory. 

Finally, we give some results of algebraic independence. What we show is that 
the elementary arithmetic functions are algebraically independent. This requires 
some deep results in algebraic geometry. 

2. The Prime Divisor Function 

Let us now study the function defined by 


(1) w(n)=^l. 

P\n 

In words, a >(n) is the number of distinct prime divisors of n. This function is 
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obviously additive. Consequently, we have the equation 



= f 2 (s)/f(2s). 


From this, we obtain the representation 
(3) 2 w < n > =£ v(k)d(n/k 2 ). 

k 2 \n 

Exercises 

1. What do we get if we write f 2 (s)/f(2s) as f(s)(f(s)/f(2s))? 

2. Show that 2 W ^ < d(n ), with equality only if n is squarefree. 

3. The Mean Value of 2 w(w) 

Using the representation above we readily find 

N 

(1) ^2 w(n) = c 1 Ar + 0(V?V). 

n =1 

We leave it to the reader to determine the constant c x . 

4. The Mean Value of c o(n) 

We have 

N N / x 

n -1 n~ 1 pin 



(i) 
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-"Si 

p<N 


+ 0 (N). 


Using the prime number theorem, we have 


( 2 ) 


V ^loglogiV + Cj +0(1). 

“ n 


p*ZN 


Hence, we have the result 


N 

(3) ^ o>(h) = AHog log N + 0(TV). 

n =1 


However, we do not need the prime number theorem. An elementary argu¬ 
ment suffices. We leave it as an exercise for the reader to supply it. 


Exercises 


1. Find the asymptotic order of 1^=1 log d(n). 

2. Find the asymptotic order of 1,^=1 log 0 (h). 

3. Show that 


oo 



4. Find the Euler product for a w ^ jn s . 

5. By differentiation, obtain 


oo 



n =1 


z. „s 

n =1 


5. The Order of Magnitude of \ o?(n) 

To evaluate this sum, we write 

N N 

2 w 2 («) = 2 w (") 

n=l n~l 



(i) 
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Inverting the order of summation and using the result of the previous exer¬ 
cises, we have 


N 

(2) 2 w2 («) = W(tog log TV) 2 + 0(TVlog logTV). 

/I = l 

Combining these results we have 
N 

(3) ]jT («(») - lo §logTV) 2 = 0(TVlog logTV). 

n=l 

From this result, we see a meaning we can attach to the statement, “Most 
numbers in the interval (1 ,iV) have log logTV prime divisors.” 


6. The Average Order of o o(g(n)) 

Let us now consider the order of magnitude 2^ <jo(g(n)) where g(n) is an 
irreducible polynomial in n . We have 

o ) 2 "<«("))“2 (2 ! ) 

n =1 n =1 p|^(n) 

To estimate this last sum, we break the second sum into two parts, 1 < JV 

and p> N. It is clear, since g(n) is a polynomial that it has only a finite number 
of prime factors p>N. Consequently, we have 

N 

( 2 ) ^ «&(»)) = 2 ( 2 0 + o(Ar) 

p<iV £<n)^0(p) 

=tv ^ ^M+o(tv), 

P<AT P 

where, as before, p(p)is the number of solutions of the congruence. 

Thus, finally, we obtain the results 
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N 

(3) «fe00) = Nlo % lo % N + W- 

n=l 

Exercise 

1. What do we get if the polynomial is reducible? 

7. The Sieve Method of Selberg 

Let us consider the following method of Selberg, which in many parts of num¬ 
ber theory replaces the sieve method of Brun. 

Suppose that we have a finite sequence of integers {a k }, k = 1,2, * • • ,TV, and 
we wish to determine an upper bound for the number of the a k which is not 
divisible by any prime p < z. 

Let {jc v }, 1 < v < z, be a sequence of real numbers such that x t = 1, while 
the other x v are arbitrary. Consider the quadratic form 

N 2 

a) q (*> = Z( 2* v ) • 

fc=l ' v|<fc ' 


Whenever a k is not divisible by a prime < z, the sum Z v \ ajc x v yields x t = 1. 
Hence, if we denote the number of a k , k = 1,2, * • • ,N, divisible by any prime 
p < z by/(TV, z) we have 

N 2 

(2) W.z)< 1(2-) 

k=\ v\a k 


for all x v with Xi = 1. Hence 


(3) 


/(TV, z) < min 

x v 


2(2 *)’■ 

k =1 v|a^ 


To evaluate the right-hand side, we write 

l 


Q(X) = ^ 

Vl,V 2 <Z 



( 4 ) 



160 


Analytic Number Theory: An Introduction 11.6 


where a k runs over the sequence a u a 2f • * * , a N . Here K denotes the greatest 
common divisor of v x and v 2 . 

Now suppose that the sequence {a k } possesses sufficient regularity properties 
so that there exists a formula of the type 

(5) y i =—+ r ( p ) 

** f(p ) 

k=l,2,-' 4 ,N 


where R(N) is a remainder term and the “density function” f(p) is multiplica¬ 
tive, i.e., 

(6) /(PiP2)=/Oi)/0°2) f or(Pi,p 2 )= 1. 

Then 


(7) 


y i = a 


nm +R (hii\ 

f(v v )f(v 2 ) \ K / 


Using this formula in (4), the result is 


( 8 ) 


m,z)< 


2 

Vl,V2<Z 


x Vl x v /(K) ^ 

f(vi)f(v 2 ) 


2 

Vl,V 2 <Z 


Consider now the problem of determining the x Vi 2 < v < z, for which the 
quadratic form 




/><*)= 2 

Vi,V 2 <Z 


f(v\)f(v 2 ) 


is a minimum. 

To do this, we introduce the function 


(10) /i («) = n(d)f(n/d). 

d\n 


If « is squarefree, we have 
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(ii) 


p\n 


Using the Mobius inversion formula, (10) yields 


(12) /(*) = £ A («)= £/i(h). 

n\K n\v i 

n\v 2 


Using this expression for f(K ) in (9), we see that 


(13) 


P(x)= ^ fi(n) 


n<z 



2 


Now perform a change of variable. Set 


(14) 




n |v 
v<z 


x v 

m' 


Then, using the Mobius inversion formula once more, we have 


(15) 


*v 

m 


-i 

n<z/v 




The problem is thus that of minimizing 


(16) 


P= ^fdnW 

n<z 


over all {.y k } such that 


(17) 


I 

n<z 


M(n>„ = 



It is easily seen that the minimizing {y k } is given by 
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(18) 


y n = 


M(«) 

h(n) 


1 


y mV) 


and that the minimum value of the form P is 


(19) 


1 


2 

n<z 


m 2 (») ’ 

/i(«) 


The corresponding x n are determined by the relations 


( 20 ) 


x„ 


. m 


y p 2 (p) 
£ifi (p) 


2 

p<z/n 


ju(p)M(p») 

/i(P") 



/«/ 


I ? 


'(/>) 


P<2 


/i(p) 


2 

p<z/n 

(P,n)=l 


m 2 (p) 

/l(p) ' 


Inserting these values in (8), we obtain the inequality 


( 21 ) 


+ ^ |x Vl x V2 *(ViV 2 /A’)|. 

y l* v 2 


If z is chosen properly, say N V2 ~ e , the upper bound is nontrivial. 

A particularly simple case is that where a k = k, so that /(p) = p. We then ob¬ 
tain a bound for the number of primes less than or equal to N. 


Exercises 


1. Prove that the number of integers less than x which have exactly k prime 
factors is bounded by 

c 1 x(loglog.x +c 2 ) k ~ x 
(fc-l)llogx 
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2. Prove that the number of integers n such that n and n + 2 are both prime, 
twin primes, less than x is bounded by c 3 x/(logx ) 2 . 

3. Hence, prove that D 1/n, where the summation is over twin primes, converges. 

4. Consider the Goldbach representation, p x + p 2 = «, where /? 1 and p 2 are 
primes. Prove that the number of such representations is bounded by c 4 w/ 
(log«) 2 . 

5. Obtain a bound for the number of primes in the arithmetic progression 
an + 6. 

6. Obtain the bound for the number of primes in the arithmetic progression 
an + b in the interval (x , x + y). 

7. Find the order of magnitude of d(p + 2). 

8. Extend Selberg’s method of the case where the quadratic form is of the 
type Z% =1 f(nX' * ) 2 . 


8. The Algebraic Independence of Arithmetic Functions 

The search for relationships between various arithmetic functions has, in the 
past, resulted in many interesting identities. These efforts have been centered on 
the problem of evaluating various types of sums of arithmetic functions in terms 
of others. In this section we will adopt a negative approach to the problem, and 
consider the question of when one can assert the algebraic independence of a set 
of arithmetic functions. 

Before proceeding to discuss what types of arithmetic functions we intend 
to consider, we shall make these terms more precise. 

DEFINITION 1. A single-valued function/(«) will be called arithmetic if it is de¬ 
fined to have real values for all positive integral values of n. 

DEFINITION 2. An arithmetic function f(n) 0 will be called multiplicative if, 
for m and n relatively prime,/(mn) =/(m)/(n). 

Perhaps the best known examples of multiplicative arithmetic functions are 
n; d(n ), the number of divisors of n; o(ri) y the sum of the divisors of n; 0(n), the 
Euler function; p(n) y the Mobius function; 2°^ where cj(n) is the number of 
distinct primes dividing n. 

DEFINITION 3. A set of arithmetic functions fi , * * * ,/)v will be called alge¬ 
braically independent over the real field if there exists no polynomial P(xj, 
' * *, x N ) =£ 0 with real coefficients, irreducible over the real field, such that 


P(fu 


> > j 




for all positive integral values of n. 
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In this section we discuss only multiplicative arithmetic functions. 

An interesting product of the investigations of this section is the following 
corollary: 

The functions n, 0(n), d(n), 2 W ^, a(n) and p( n) are algebraically independent. 

The proof of this result is circuitous in that it seems essential to discuss alge¬ 
braic relationships between general multiplicative arithmetic functions first. 

Since the polynomial relations between two multiplicative functions can be 
discussed more fully than those involving more than two functions, this case 
will be discussed first. The extension of some of the results to more than two 
functions depends upon the results obtained for the case of two. 

To begin with, it should be noted that any two arithmetic functions /(/?), 
g(n ), each of which takes on only a finite number of different values are alge¬ 
braically dependent. This follows from the fact that there are only a finite num¬ 
ber of points [/*(«), g(n)] in the plane, and by Bertini’s theorem an irreducible 
algebraic curve can be passed through the points. 

DEFINITION 4. A function f(n) which assumes only a finite number ot differ- 
erent values for all positive integral n> will be said to be finitely valued. 

DEFINITION 5. A function f{ri) which assumes an infinite number of distinct 
values on a set of numbers of the form p a *, for only a finite number of primes 
p, will be called singular. 

DEFINITION 6. A set of functions /*(«), 2 < k < N, will be said to be singular 
if all the f k (n) simultaneously take on an infinite number of distinct values only 
on a set p a fc, where p is a fixed prime for all Jc. 

Also to be noted at this point is the fact that if two arithmetic functions fin) 
and g(n) constitute a singular set, they may be algebraically dependent, even 
under the added condition that / and g be multiplicative. To show this we make 
use of the fact that a multiplicative function may be defined arbitrarily at the 
values n = p k , fc=l,2,-**,pa prime. Consider any irreducible algebraic curve 
through the points [0, 0] and [1,1] and for some fixed prime p let the couples 
[/(p a 0> g(p a 01 correspond to an infinite set of points on the curve. For all 
other n ¥= 1 define f(n) and g(n) to be zero. These functions are then clearly 
both multiplicative and algebraically dependent. 

We can moreover give an example of three multiplicative arithmetic functions, 
forming a singular set, which are positive, and still algebraically dependent. 
Suppose that [a if b it c t ] are the finite number of sets of values assumed by 
[/> g, h] for the values n =£ p a *, (n, p) = 1. Then let all the points [ f(p a i ), g(p a '), 
hip “*)] lie on an algebraic curve C through the points [a if bc t ] and consider 
the linear transformations 


a t 0 0 

tf 0 b t 0 

_0 0 Cj_ 
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Operating with the inverses of these we get a set of curves C% = t]\C ). Then by 
Bertini’s theorem we can pass an irreducible algebraic surface S through the 
curves C and C,*. The surfaces tj(S) will all intersect in the curve C. Since the 
points [/, g, h] he on S, f g, and h are algebraicahy dependent, and the above 
construction insures complete consistency with the multiplicativity of these 
functions. 

The above remarks indicate the nature of some of those cases of multiplica¬ 
tive functions which are algebraically dependent. They also serve to establish 
somewhat the necessity of the hypotheses which will be found in the following 
theorems. 

We will now state some auxiliary results which will be needed later. 

LEMMA 1. If f(n) is not finitely valued there is no polynomial relation P[f(n)] 
= 0 . 

This is clear. 

LEMMA 2. If f(n) is a multiplicative arithmetic function , then either (a) f(n) 
is finitely valued , or (b) f(n) is singular , or (c) there exists a sequence of integers 
5 k relatively prime by pairs such that f(5 k ) takes on an infinite number of dis¬ 
tinct values. 

The proof of this follows easily by showing that the negation of (a) and (c) 
implies (b). 

THEOREM 1. If f(n) and g(n) are two multiplicative arithmetic functions such 
that [f, g] is not singular ; and at least one of the two is not finitely valued , then 
f and g are algebraically independent , unless f = g r where r is a rational number. 

PROOF. Suppose that / and g are algebraically dependent, so that P[f g] = 0 
where P(x, y ) is an irreducible polynomial. 

If f(n) is finitely valued, then g(n) takes on an infinite number of different 
values on a set n h i = 1, * • •. We write 

p[f,g] =p m [f]g m +-+p 0 [f]- 

Then since / is finitely valued there is an infinite subsequence of the n t such that 
f (rif) = A, a constant, and £'(/!,*) assumes an infinite number of different values on 
this subsequence. But then 

P m {A)tr+---+P 0 \A] =0 

would have an infinite number of roots, which in turn implies that 
P k [A] =0, (A; = 0, * • •, m). 


Hence 


P k (x) = (x-A)Q k (x) 
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and 


P(x, y) = (x-A)Q(x, y), 

contrary to the assumption that P(x, y ) is irreducible. 

On the other hand if both / and g are not finitely valued, Lemma 2 implies 
that there exists either (a) two infinite sequences {a,-}, {ft} such that (a*, ft) 
= 1 for all i and /, and at least one of the two functions, say g, takes on an in¬ 
finite number of distinct values on both {a*} and {ft}; or (b) two primes p and 
q such that p ± q, and / takes on an infinite number of distinct values on a se¬ 
quence p ai 9 and g takes on an infinite number of distinct values on a sequence 

Now we have P[f(rt), #(«)] = 0, and if (m, ri) = 1, f(m) = a, g(m) = b, 

P[af(n), bg(n )] = 0. 

Case (a). We fix the couple (a, b) as any element of the sequence [/(a,), 
£■(<*;)] and let n run through the sequence {ft}. Then P[f, g] = 0 and P[af, bg] = 
0 have an infinite number of roots in common. Eliminate / between these two 
equations thus obtaining the resultant R(g ), which must vanish for the infinity 
of values #(ft). This clearly implies that R(g) = 0, since R is a polynomial in g. 
Since P[f, g] is irreducible we must have 

P[af, bg] = c(a, b)P[f, g], 


or 


^Ai^bifg’ = c(a, b^Aijfgi, 

and thus a l b J = c(a, b). Since P[f f g] is irreducible and must contain both vari¬ 
ables it must have more than one term. Hence, we get an equation of the form 
a l b ] = a k b f which we can write as 


a r b s -1=0 or a r -b s = 0. 

Thus in either case Q(x, y) = x r y s - 1, or x r -y s 9 has a degree less than or equal 
to that of P(x , y). Now since (a, b) can run through the infinity of couples 
[/(«/), £(<**•)], Q(x, y) and P(x, y) have an infinite number of roots in common. 
Eliminating x, we again get that the resultant /?!(» = 0. But this implies that 
P(x, y) divides Q(x, y), and since the degree of Q(x, y) is less than or equal to 
the degree of P(x, y) we must have P(x, y) = cQ(x, y ), and hence Q(f\ g) = 0 
which is the conclusion of the theorem. 
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Case (b). The argument is much the same as in Case (a), (a, b ) is first fixed 
in the sequence [fip 01 *), £Cp a 0] and we let n run through the q&. Then as be¬ 
fore we get Q(x, y) to have the same forms as noted above, having the roots 
[f(p«i), g(p a *)] in common with P(x, y). Then proceeding from this point we 
now eliminate y between Q(x, y) and P(x, y) obtaining the resultant R 2 (x). 
But since R 2 [fip* 0] = 0, and f(p ai ) takes on infinitely many values, R 2 (x) = 
0, and the conclusion follows as before. 

COROLLARY. If f(n) is a multiplicative arithmetic function , and n and f(n) are 
algebraically dependent , then f(n) = n r for some fixed rational number r. 

Before proceeding to the general case of any number of multiplicative arith¬ 
metic functions, there is one more pathological case to mention produced 
by our definitions of r. For N > 2, we can have a set of arithmetic functions 
such as the following: 


/l(Pl) = C ai 

*0, 


/i(n) = 0 

otherwise; 

/l(P2> = C<x 2 

*0, 

/ 2 ( 1) = 1, 

> 

Is" 

'W' 

II 

o 

otherwise; 

/a Of) =Coc 3 

*0, 


f 3 (n) = 0 

otherwise. 


If the Pi t p 2 , P 3 are different primes and the c a . are any infinite sequence of 
real numbers, f x , f 2 and / 3 clearly constitute a nonsingular set of multiplicative 
arithmetic functions. Also, they are algebraically dependent since 

Q\fl,h.h\ =/l/ 2 +/2/3-2/ 1 / 3 , 

is an irreducible form, and vanishes for all n. 

For obvious reasons our hypotheses will have to be formulated so as to elimi¬ 
nate such sets of functions from consideration. Due to the variety of peculiar 
combinations of functions which can arise, as is indicated by the above discus¬ 
sion, we will not be able to prove as sharp a result for the case N > 2 as in 
Theorem 1. However, we can obtain a workable result. 

Theorem 2. Let f k (n), k = 1,2, * • •, N, be a set of multiplicative arithmetic 
functions. If a polynomial P(xj, • • •, x N ) ^ 0, reducible or irreducible , vanishes 
on the sets (fi(n), • • • , f N (n), n = 1, 2, • • • , then for every sequence of pair¬ 
wise relatively prime integers {nj there exists a subsequence {nj} on which there 
holds a relation of the form 

(i) ft 1 ■ ■ -f N * N =/f‘ • ■ -f N pN , 

where the cq, ft are nonnegative integers and , for at least one i, ^ft. 

In (1) it is understood that for k = 0, f k = 1. Also a relation 
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(2) A* 1 ••■^ = 0 

on {nj} is of the form (IJ, i.e., ffi • • -fgN ... fg*N {nj}. 

PROOF. Suppose that P(/i(«), • • * , f N (n)) = 0 for all Let {«,} be any se¬ 
quence of relatively prime integers. If P has only one term, the theorem is 
proved by the note above. If P has more than one term, write P(fi(n ), • * * , 
Ar(«)) = Af* 1 (n) • • • f]tf N + Bff 1 • • • + ♦ • •. Either (1) holds on {wjor 

there is an m e {«/} such that a = /f 1 (m) • • • j f£ N (m) ^f/ 1 (m) • • • j= b. 
In the second case, P(/i(mw), • • • J N (mn)) = Aaf* 1 (n) • • -jj^Cii) 

* 0 *JSv iv («) + * ■ •andoP(/ 1 (/i), • • *,/Ar(«))-P(/i(^«), " • ,f N (mn)) = B(a-b) 
fi 1 («) *** fJj N (ri) + ■ * * is not identically zero, a ^b, vanishes for all n such 
that (n, m) = 1, and has fewer terms than P. It is clear that either (1) holds on 
some {n'i } or the process can be continued until the resulting polynomial is of 
the form (2) and vanishes for all n with (n, mj = 1, i = 1, • • • , k, k< the 
number of terms of P. This completes the proof. 

We will now apply the above results to certain rather well-known sets of arith¬ 
metic functions. 

Theorem 3. If f k+1 = S d|n f k (d), k = 0, ♦ • • , N, f 0 (n) is multiplicative and 
f 0 (p) takes on an infinite set of values over some set of primes, then f 0 , * • • , 
f N are algebraically independent. 

PROOF. Since f 0 (n) is multiplicative, all th ef k (n) are multiplicative. Also since 
fo(p) takes on an infinite number of distinct values on the sequence of primes, 
so does fi(p) = fi-\(p) + 1 =fo(p) + *. Thus if/ 0 , * ■ * ,/at were to be algebraically 
related we would have 


[/oCrt + ii]* 1 ■••[/o(P) + y^=[/o(P) + ^ + i] afc+1 •••[/o(P) + k] 0 *, 

where ij =£ ij , 0 < < N, a,y > 0. But since / 0 (p) takes on an infinite num¬ 

ber of distinct values this would require that we have 

(3) [x + ii]* 1 "[x + i k ] < *K = [x+i k+1 ] 0l K+ 1 ••• [x + l R ] a R. 

But this is impossible, for one side has the root x = -i t , and the other side does 
not. Thus by Theorem 2, / 0 , * * * ,f N are algebraically independent. 
COROLLARY, n, 0(n), a(n) are algebraically independent. 

PROOF. Take/ 0 = <p(n) in the theorem. Then f x = n, f 2 = o(n). 

THEOREM 4. n, 0(n), a(n), d(n), 2 ^ are algebraically independent. 
PROOF. We first show that there exists an infinite sequence of primes satisfy¬ 
ing the following five conditions: 

(1) Pi = -1 (mod 3), pi> 3; 

(2) pi £ 1 (mod p f ); 

(3) p,*2 a + l; 

(4) p\ + Pi + 1 £ 0 (mod pf) and mod 3); 
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(5) pj = 1 (mod 2), Pj- 1=0 (mod/?/) implies that pf + p t + 1 £ 0 (mod/?/). 
We begin by choosing p x as any prime of the form 3n - 1 which is not equal to 
2 a + 1. Suppose then that p x , • * • , p r have already been constructed. We will 
now establish the existence of a suitable p r+1 . Let 2 6 * be the highest power of 2 
dividing /?,-!, and set 


7 = 3fJ(p, ? +Pi + l)Pi{(Pi - l)/2 6 '}. 

r-\ 


If we consider the form yx - 1, Dirichlet’s theorem gives us a prime of this 
form, which is greater than all the p t and all the prime divisors of pf + /?/ + 1, 
i* = 1, • • • , r. We let this prime be /? r+1 . Clearly p r+1 satisfies (1). That it can be 
chosen to satisfy (3) also follows easily from the fact that E/?' 1 • log p (p = 
-l(mod7)) diverges whereas 'Lp ml \o%p{p = 2 a + 1) converges. (2) follows 
since 


Pr+i -1=7*- 2 =-2^0 (mod/?/), 

where i = 1, • • ,r. It is obvious that pf +/?,* + 1 # 0 (mod p r + x ) and pf + x + /? r+1 
+ 1 = 1^0 (mod 3) so that (4) is satisfied. Finally if sy is any odd divisor of 
Pj- 1, )<N, 


p r+i = -1 (modsj) 


so that /? r+ x - 1 = -2 ^ 0 (mod sy). 

P?+i + Pr+i +1 = 1 £ 0(mods/). 

In addition if s r+1 = l(mod2), s r+1 |/? r+1 - 1 and s r+1 \pf + /?,- + 1, these 
imply s r +117 and hence /? r+1 =-1 = 1 (mods r+1 ), which is impossible. This 
gives us the desired infinite sequence of primes Pi, * * * , p r , Pr+i , * * ‘ • 

We now form from these primes the numbers 


S*-Pl, 52“P2P3 j 


2„2 


«*=P ?+1 •• -Pr+k, r = 


_ fc(fc - 1) 


Since the functions simultaneously take on an infinite number of distinct values 
on the sequence {§£}, and since the 8* are relatively prime to one another, we 
are now in a position to apply Theorem 2. Suppose then that 


«“10(rt)f 2 (o(«»“ 3 (rf(««))“ 4 2 w(n) “ s = (<p(n)f 2 (o(n)fi (d(nn)f* 2 0J <"^ 
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over some subsequence of the { 8 k }. This gives 

(np 2 )“i(np(p- i))“ 2 (n(p 2 +p +1 ))f 3 3 ta * 2 tas 
= (np 2 /Knp(p- i)/ 2 (n(p 2 +p + i)p3 k04 2 k0s . 

From the manner in which the p’s were chosen, each p t - 1 contains an odd 
prime divisor greater than 3 which divides none of the other terms which appear. 
Hence a 2 = & 2 . Then since p t > 3, and divides no term of the form pf +pj + 1, 
we must have a x = 18 !. Similarly since pf + p t + 1 # 0 (mod 2), (mod 3), we get 
of 3 = j 83 . This leaves 3 0£4 2 Q!5 = 3^ 4 2 ^ 5 so that a 4 = p 4 and a 5 = j 3 5 . Thus by 
Theorem 2, the set of functions under consideration are algebraically independent. 

Corollary. The set of functions n, 0(n), a(n), d(n), 2 w(n) , ju(n), are algebra¬ 
ically independent . 

[p(«) is the Mobius function. Since it has a range of only the three values 0, 
- 1 , 1 , its algebraic independence of other sets of functions cannot be obtained 
directly from Theorem 2.] 

PROOF. Let P{f x , • • • , / 6 ] = 0 be an irreducible polynomial relationship be¬ 
tween the functions mentioned in the corollary. If we let / 6 = p(n) we may write 

P* [ti*)YPr[fu- • ’ ,/S + ’ • ‘ + M(«)Pl \fu * * * ./*] +^0[/l, * * * ,/ 5 ] = 0. 

But over the sequence {§£} constructed in the proof of the theorem = 0. 
Thus over the sequence constructed in the proof of the theorem we must 
have P 0 [f x , • • • , f s ] =0. But this in turn is impossible by Theorem 4. Hence 
P 0 [fu"-,fs] = 0 which is a contradiction since then P[f x , • * • ,/ 6 ] is divisi¬ 
ble by f 6 and hence not irreducible. 

Exercise 

1. Prove that the assumption that f(n) is multiplicative is unnecessary. 
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12. Tauberian Theorems 


1. Introduction 

In this chapter we wish to consider some Tauberian theorems. Tauberian 
theorems may be thought of as a method for converting from one average to 
another, and thus are an indispensable tool for the analytic number theorist. 

In the first few sections, we give the fundamental results of Hardy and Little- 
wood. Then we apply the continuous version of this result to the error term for 
0(h) and the mean value of £ 2 (s) along the critical line. 

In the final section we say a few words about the Erdos-Selberg method in 
the theory of primes. 

2. The Tauberian Theorem of Hardy and Littlewood 

It is elementary to go from the asymptotic order of a n to the behavior 
of a n x n as x -> 1. A result of this type is called Abelian. To go in the 
reverse direction we require some condition on a n . Hardy and Littlewood ob¬ 
tained many results of this type. Basically, we can pass from the behavior of 
a n x n as x 1 to the behavior of %n=i a n if we know that a n > 0. 

A fundamental result of Hardy and Littlewood is that if 


( 1 ) 


n=l 



as* -► 1 
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and a n > 0, then we have 


( 2 ) 

Exercises 



^ 7V°(logA0 & 

r(« + i) ’ 


a> 0. 


1. Prove that 

©o oo 

(1 -x) J a n x n = ^(a„- a „.i )x n , a-j = 0. 

n=0 n=0 

2. Hence show that we can use a Tauberian theorem if we add the hypothesis 
that a n is monotone increasing. 

3. Write 

oo 

What is the significance of the coefficient p„? 

4. Shew that p„ is monotone increasing and use an appropriate Tauberian 
theorem to obtain the asymptotic behavior of p„. 

5. On the basis of the prime number theorem, find the behavior of S p x p as 
*-► 1 . 

6. Show that 


^(x p f b n x n , 

p n~l 

where b n is the Goldbach number, the number of solutions of p v + p 2 " n > 
where p x and p 2 are primes. 

7. Find the order of magnitude of E^ =1 b n as N-+ °°. 

3. Continuous Version 

The continuous version of the foregoing result is also useful. If we have 

Jo y 


(i) 


asy-^0. 
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and /(jc)>0,then 

(2) [ T f(x)dx~ Ta(logr)b a > 0. 

W Jo r(a + l) 

Exercises 

1. Consider the equation u(x) = f(x) + Jq k(x - y)u{y)dy, commonly called 
the renewal equation. Show that u(x) is nonnegative if f(x) and k{x) are non¬ 
negative. 

2. Show that this equation can be solved in terms of the Laplace transform of 
u(x). 

3. Obtain conditions on f(x) and k(x) which permit the determination of the 
asymptotic behavior of u(x)dx. (This equation plays a very important role 
in the theory of branching processes. See the book 

Mode, C. J., Multitype Branching Processes , American Elsevier, New York, 
1971. 

Other information concerning the renewal equation can be found in 

Bellman, R. and K. L. Cooke, Differential-Difference Equations , Academic 
Press, New York, 1963. 

4. The Mean Value of (1/2 + it) \ 2 

We begin with the fundamental representation 

(1) r(s)f (s) = 

Jo e x -l 

Regarding this as a Mellin transform, we have the inversion formula 

(2) -j— =-L [ x-*r(s)Hs)ds. 

e x - 1 2m Jc 

Here C is a straight line parallel to the imaginary axis to the right of the line 
o- 1. 

Let us now shift the contour. If we shift past the line o = 1, we encounter a 
simple pole at s = 1. This result is the formula 


( 3 ) 


- -=~ f x-°r (s)r(s)ds. 

e -1 x 2m Jc x 
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The contour C\ is a line a = It remains to choose x adroitly. Let us choose 
the value x = ie l °. If we use the fact that i = e^ 2 we see that (3) is a Fourier 
transform. 

We now use the Parseval-Plancheral formula. The details are given in a mono¬ 
graph by Titchmarsh cited at the end of the chapter. 

5. The Error Term in d(n ) 

Let d(n) denote the number of divisors of n , and consider 


( 1 ) 


D(x) = ^ d(n). 

n<x 


It was proved by Dirichlet that 

(2) D(x) = x log* + (2C - l)x + A(x) (x -+ ~), 

where A(x) = 0 (x^ 2 ), and C is Euler’s constant. Subsequently, it was shown 
by Voronoi that A(x) = CHx^logx), and the estimate has been continually im¬ 
proved since, although the precise result is still unknown. 

In the other direction, it was proved by Hardy that for some constant k y 

(3) |A(x) | > kxV* 
for an infinity of values of x -> and even that 

(4) |A(jc)| > k(x log jc) 1 / 4 log log* 

for an infinity of x -► °°. A result of this type is called an f2-result, and (4) is 
written 

(5) A(x) = Sl(x log x) 1/4 (log log x). 

If we are interested, not in exceptional values, but in the average value, as in 
the case of d(ri) itself, we have the following result, also due to Hardy 

(6) J^\A(x)\dx = 0(T /i+e ) (T-k> o). 

The result we wish to prove in this chapter is a refinement of this result, 
namely 
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(7) r^fdx-CilogT (r->~), 

J 1 X* 

where c x ^ 0 is a constant. Hence A(x) = ). 

The same method as used in the proof of this theorem yields a similar result 
for the error term in the Gauss circle problem, which concerns itself with the 
asymptotic order of 


( 8 ) 


^ r(n) = R(x) 
n<x 


where r(n) is the number of solutions in positive integers of a 2 + b 2 = n. 
We shall consider the error term in the estimate of 

(9) 

n<x 

where 0(«) is the Euler 0-function, and derive a result similar to that above. 

We mention in passing that this method yields the known results for the func¬ 
tions tt(x) - Li(x ) and \p(x) - x as far as 12-results are concerned, but not the 
mean-value theorems, which seem to require more knowledge concerning the 
distribution of the zeros of the Riemann f-function (see 

Ingham, A. E., The Distribution of Prime Numbers , Cambridge University 
Press, London, 1932. 

The method used in this section for dealing with problems of this type seems 
due to Titchmarsh. 

We have 


(10) ? 2 ( s )=2^r (**(*)>!)• 

n=l 

Using Perron’s sum formula, we deduce 

(11) y, ^( M ) = — I f 2 (w) — dw (x > 0, nonintegral). 

n< x 

Since the integer values of x form a denumerable set, we can ignore the be¬ 
havior of the sum-function at these points. 
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We now shift the line of integration to the left of the line 2 + it until we reach 
the line a + it, X A < a < H. This is accomplished by the standard method of first 
considering the rectangle whose vertices are 2 ± iT, a ± iT, and then allowing 
T to tend to *>. There is no difficulty in justifying this procedure using the 
order relation |f(% + zY)l = 0(\t \ % ) as t -► °°. The only singularity of the inte¬ 
grand is at w = 1, and the residue at this point 

(12) R(x) = x\ogx + (2C-l)x. 

Thus 


( i3) ^ m . Kx)= ±.rqs^.^ dt . 

^ 2 m J (a + it) 


n<x 


. Let us now make the substitution x = e u , so that the right side is a Fourier trans¬ 
form. Thus 


(14) 


*A(e“)=— f°° ^ (a + it) e iut dt. 
2 in J-oo a + it 


Before proceeding further, we shall show that f 2 {a + ifjfa + it)" 1 belongs to 
Z, 2 (-°o, °°). Once this has been established we can use the Parseval-Plancherel 

theorem to obtain the following equation: 


(15) 


f°° e~ 2au A 2 (e u )du = — f ~ ^ (a + lt ^ dt 
J-oo y 2 t rJ-oo a 2 +1 2 


= I f°° l£(g + ft)r rff 
IT J o a 2 + 1 2 

To prove the L 2 -result, we require the following three results. 

( 16 > J*\s(p+it)pdt~c 2 (pyr, i >o>Vl (r+~). 


(17) 


f(s) = 2 V _ 1 sinyr(l -s)f(l -s). 


(18) |r(l -a-if)|~C3(fl)|/| v 4 " a e " w|fl / 2 (|f|-*~). 


Combining (18) and (17), we obtain 
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(19) 


|f(a + /Y)|~£ , 4 (flr)|?| % '°li'(l -a-it )| (|/|-»•«>). 


Thus the convergence of the right side of (15) is dependent upon the con¬ 
vergence of 


( 20 ) 


/; 


f 4( *- a) |f(l - a-it)I 4 


a 2 + f 2 


dt. 


which in turn is convergent if 


( 21 ) 


f” ir(l-a- 

Ji f 40 


-«0I 4 




We now integrate by parts, transforming (21) into 


( 22 ) 


/>>-<■- 


Jki 


it)\ 4 dt~r 

—w; 




i'0l 4 df 


,l+4a 


dt. 


Using (16) we see that the integrated part vanishes, and the integral is less than 


(23) 


Cl 


J i r 


tdt _ c 3 (a) 


1+4fl a — 1/4 * 


provided 0 > Va and < Vi. This permits us to assert that (15) is valid. 

We now retrace our steps in the preceding proof, and establish the asymptotic 
relation 


P4) /-!(^*). 

Jo t a- 1/4 

Once this has been established, we shall have the following relation 

(25) f“ e * 2 ( a -i / 4 )“ e '“/2 a 2 (e“) - (a -> %). 

J.°o a-1/4 

It is easy to see that the integral over the interval (0, -°°) is finite as a -► 
and thus 

Jo 


( 26 ) 
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or 


(27) | e' m e u > 2 A 2 (e u )du ~ c 7 /s (s -*■ 0). 

Jo 

The use of a familiar Tauberian theorem, due to Hardy and Littlewood, (27) 
leads to 

(28) ( T e' u ' 2 A 2 (e u )du ~ c n T (T + ~), 

Jo 

whence a change of variable yields (7). 

We now have to prove the relationship given in (24). It is again sufficient to 
consider the behavior of (21). It is now important to note that since (16) is 
valid for % < o < 1, we have 

(29) J T If (1 - a - it)\ 4 dt = c 2 (J)r + eT (e small), 

for a close to % and T large. Here we have used the fact that c 2 (a) depends con¬ 
tinuously upon a, which follows from the form of c 2 (a). 

Using (29) and integrating by parts as in (22), we obtain the desired result 
(24) and the proof is complete. 

We now consider the function defined in (9). The ^-result we obtain here, and 
more, was obtained by Chowla and Pillai, using very different methods, but the 
mean-value theorem seems to be new. 

Consider the generating function for the Euler 0-function, 


(30) 


K«-i ) -f#) 

f(j) " k s 

fc=i 


(R(s)> 2). 


Applying the Perron sum-formula, 


(31) 


2 <K«) 

n<x 


J_ f 3+i °° f (S - 1) ds 
2m J 3 -/« f(s) s 


Shifting the line of integration to the line 1 + a + it, 0 < a < 1, we obtain the 
equation 
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(32) 


<Kn) 


3x 2 


n< x 


2 7T 



S(a + it) x 1+aHt dt 
f (1 + a + it) (g + it) 


Denoting the left side of (32) by Aj (x), we obtain as before 


(33) 


f °° e - 2<fl+l »A\{e u )du 


J_ f°° jjg+ it) 2 dt 
2n J.=o f(a + 1 + if) a 2 +t 2 ' 


Using (19) and a bound on If (1 + a + it)!' 1 , we have 


(34) 


Ka + ft) 2 

f(l +a + it) 


= o(t 1 - 2a 0ogki 4 )). 


Thus the right side of (33) exists and, using Plancherel’s theorem, (33) is valid. 

To determine the asymptotic behavior of the left side of (33) as a -> 0, we have 
to consider the behavior of 


(35) 



If(l-g-ft)! 2 dt 
|f(l +a + it)\ 2 t 2 ** 1 ' 


We expect that as a -» 0, this acts like 


<36) /TA-? ( ^ 0) - 

To show this, we use the following equation. 

/; 

= r ~ + r (i^o-a+tpi 2 -if(i +g+ft)i a ) _ dt 

J i t 2a+1 J i |f(l + g + if)| 2 t 2 ^ 1 


(37) 


|f(l -g + /f)| 2 
|f(l +a +i'f)| 2 t 20+1 


-£♦/. 

•[/; 


l+a 

1-a 


{f (1 - a - »QfV + it) + f (1 + a + /Qf'(a - ft )} 
t 2 fl+ 1 lf(l + g + ff)l 2 




Consider the integral; it is less than 
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(38) 



iMi-«-ft)iir(g+ft)i A i 
t 2a+1 1?(1 + a + it)\ 2 J 


da. 


which in turn is less than 


(39) 



f(l -a-it) 2 dt + f°° 

f(l +a + it) 7 s * J l 


f( 1 +a + it) 


r ]da. 


Applying Schwarz’s inequality again, we obtain the following bound for the 
integral: 


(40) 


r [/; 


|f(l -a-it)\ A 

tla+1 


dt + 2 


/: 


dt 


|f(l + a + it)\ 4 t 2a+1 


/*« 

J 1 


*2<z+l 


dt do. 


It is easy to show that the mean values 

UmL ( T \S(l-a-it)\ 4 dt, 
T-+<*> T J \ 


(41) 


lim — " 


dt 


t-+'~T i lf(l+fl + ft)l 4 ’ 

Urn [ T \?(o + it)\ 4 dt 

T-+°o T J \ 


exist, uniformly in a and o for a positive and close to 0 and a close to 1, and 
thus that the inner integrals of (40) are uniformly bounded by cfa as a 0, 
where c is a positive constant. Thus the entire expression of (40) is bounded by 
a constant as a 0. This proves the relation stated in (36). 

Thus we obtain the following result: 

Let 


(42) A 2 (x) = ^ 3x 2 /n 2 . 

n<x 

Then 


f T *2&ldx~clogT (r-►«>). 

J1 X 


( 43 ) 
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Hence, A 2 (x) = £2(x) as x -> °°. 

Chowla and Pillai obtained the following results. 

(44) J T A 2 (x)dx ~cT 2 
and 

(45) A 2 (xr) = £2(x log log log*) (x -* °°). 

In conclusion, we may remark that the method illustrated above furnishes a 
systematic means of obtaining both asymptotic results for the average of the 
error term and 12-results. 

6. The Mean Value of |f (H + #)| 4 

The procedure used in Sec. 4 may be used to obtain the mean value of 
If (Vi + z/)| 4 naturally, the analysis is more formidable. 

As we point out in the references, the mean value has been obtained by several 
authors using different methods. 

For the method followed in Sec. 4, we require a functional equation due to 
Wigert. References to Wigert’s original paper and a much shorter proof by Lan¬ 
dau are given at the end of the chapter. 

7. The Erdos-Selberg Method for the Prime Number Theorem 

Let us say a few words about the famous result of Erdos-Selberg concerning 
the prime number theorem. Many mathematicians have believed that it was 
impossible to find a proof of the prime number theorem without leaving the 
real line. 

The Erdos-Selberg proof consists of two steps. First, classical methods are used 
to derive a nonlinear recurrence relation. Then a Tauberian theorem is used to 
obtain from this relation the prime number theorem. 
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Miscellaneous Exercises 

1. By using the relation f R f (s)f (1 - s)ds = 0, where R is a suitably chosen rec¬ 
tangle, can we derive the asymptotic behavior of Jq If (Vl + it)\ 2 dtl 

2 . Consider the series 

V (d(n) - log w - (2c - 1)) 
n=1 

so that this series converges in the critical strip and obtain an estimate for the 
line of convergence. 

3. Show that this can be used to obtain an estimate for f (M. + it). 

4. Can the higher divisor functions be used in the same way? 
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ANALYTIC NUMBER THEORY: An Introduction 

The purpose of this book is to provide an introduction to analytic number theory. It 
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functions of a complex variable and can be used for a one-semester course, or self-reading. 




